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The real part of e' /e measures direct CP violation in the decays of the neutral kaons in two 
pions. It is a fundamental quantity which has justly attracted a great deal of theoretical as 
well as experimental work. Its determination may answer the question of whether CP violation 
is present only in the mass matrix of neutral kaons (the superweak scenario) or also at work 
directly in the decays. After a brief historical summary, we discuss the present and expected 
experimental sensitivities. In the light of these, we come to the problem of estimating e'/e in 
^ ^, the standard model. We review the present (circa 1998) status of the theoretical predictions 

^ ■ of e' /e. The short-distance part of the computation is now known to the next-to-leading order 

' in QCD and QED and therefore well under control. On the other hand, the evaluation of the 

hadronic matrix elements of the relevant operators is where most of the theoretical uncertainty still 
OO resides. We analyze the results of the currently most developed calculations. The values of the Bi 

^N) , parameters in the various approaches are discussed, together with the allowed range of the relevant 

combination of the Cabibbo-Kobayashi-Maskawa entries Im VtdVf*^. We conclude by summarizing 
' and comparing all up-to-date predictions of e'/e. Because of the intrinsic uncertainties of the 

^ ' long-distance computations, values ranging from 10~* to a few times 10"'' can be accounted for 

' in the standard model. Since this range covers most of the present experimental uncertainty, it 

, is unlikely that new physics effects can be disentangled from the standard model prediction. For 

. updates on the review and additional material see http://www.he.sissa.it/review/. 
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I. WHAT e'/e IS AND WHY IT IS IMPORTANT TO KNOW ITS VALUE 



A CP transformation consists in a parity (P) flip followed by charge conjugation (C). It was promoted (Landau, 

1957) to a symmetry of nature after parity was shown to be maximally violated in weak interactions (Wu et al., 1957). 

Until 1963 the CP symmetry was thought to be exactly conserved in all physical processes. That year, J.M. 
Christenson, J.W. Cronin, V.L. Fitch and R. Turlay (1964) announced the surprising result that the CP symmetry 
was indeed violated in hadronic decays of the neutral kaons. 

In order to interpret the experimental evidence we must consider the strong Hamiltonian eigenstates and its 
CP conjugate as an admixture of the physical short-lived Ks component — which decays predominantly into two 
pions — and the physical long-lived component — which decays predominantly semileptonically and into three pions. 

The two and three pion final states are, respectively, even and odd under a CP transformation. Therefore, in the 
absence of CP violating interactions, we would expect the Ks,l mass eigenstates to coincide with the states 

= {K° + K'')/V2 

K2 = iK° - ^0)/\/2 , (1.1) 

which exhibit a definite CP parity, even and odd respectively (we choose CP \K°) = 1-^°))- 

What was observed in 1963 was that also Kl decays a few times in a thousand into a two-pion final state, and 

accordingly that the CP symmetry is not exact. 

The violation of CP in Ks^l decays can proceed indirectly, via a mismatch between the CP eigenstates ^ ^'^'^ 

the weak mass eigenstates Ks,l introduced by a CP- violating impurity in the K^-K° mixing, and/or directly in the 

decays of the CP eigenstates. Both effects are usually parameterized in terms of the ratios (for a recent theoretical 

review see, for instance, (de Rafael, 1994)) 

^00 = / nni r^N (1-2) 



and 



{7r°Tr°\Cw\Ks) 



- {■.+,^-\Cw\Ks) ' ^'-"^ 
where Cw represents the AS = 1 weak lagrangian. Eqs. (1.2) and (1.3) can be written as 

jyoo = £ 1= ^e-2£ , 

1 — wv2 

where the complex parameters e and e' parameterize indirect (via K1-K2 mixing) and direct (in the Ki and K2 
decays) CP violation respectively. The Ks,l eigenstates are given by 



Ks = 



= ^|±i£i , (1.5) 

where £ is a (complex) parameter of order 10""^ which depends on the chosen CP phase convention. The Ki — K2 
mixing parameter e is simply related to the observable parameter e in eq. (1.4) (see eq. (1.15) below). The parameter 
ui measures the ratio: 



{{TnT)(i^2)\C-w\Ks) 



((7r7r)(7=o)l'^^w|^s) 



1/22.2, (1.6) 



where 7=1 and 2 stand for the isospin states of the final pions. For notational convenience, we identify in the 
following uj with its absolute value. The smallness of the experimental value of w given by (1.6) is known as the 
A7 =1/2 selection rule of — > tttt decays (Gell-Mann and Pais, 1954). 
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In terms of the Ks,l decay amplitudes, the CP violating parameters e and e' are given by 



and 



{{-k-k)(i=o)\Cw\Ks) ' 



£ \ {{nn)l=2\Cw\KL) {{^^)l=2\Cw\Ks) \ 



(1.7) 



(l.E 



From cqs. (1.5) (1.8) one sees that direct CP violation arises due to the relative misalignment of the Kg and Kl 
/ = 0, 2 amplitudes and it is suppressed by the A/ =1/2 selection rule. 

According to the Watson theorem (Watson, 1952), we can write the generic amplitudes for and to decay 
into two pions as 

{{'^'^)(i)\^w\K'-') = -iAi exp {i6i) 

((7r7r)(7)|£i^|^°) = -iA*j exp {i6i) , (1.9) 
where the phases Si arise from the pion final-state interactions (FSI). Using eq. (1.9) and the approximations 

|£|Im^o < ReAo , < 1 (1-10) 

the e' parameter in eq. (1.8) can be written as 



where the parameter u) can be written as 



^ = — IT ■ (1-12) 
Re An 



By decomposing the AS" = 2 weak lagrangian for the K^-K^ system in a dispersive and an absorptive components 
as M — i r/2, where M and F are 2x2 hermitian matrices {CPT symmetry is assumed), one obtains for e the 
expression 

. ^ if> f ImMi2 Imyln \ 



V AMls RcAo 

where AM^s is the mass difference of the K^-Ks mass eigenstates, M12 is the Ki — K2 entry in the mass matrix, 
and 

0e = tan-i (2AMis/ArsL) ^ . (1.14) 

In obtaining eq. (1.13), in addition to the approximations of eq. (1.10), the experimental observations that AMls — 
F5/2 and Tl <^Ts have been used. With the above approximations one also obtains a simple relation between the 
observable parameter e and the phase-convention dependent parameter e, 

Im An 

e = e + i—^. (1.15 
Re Ao ^ ' 

For detailed discussions on the role and implications of the phase conventions we refer the reader to the reviews of 
(Chau, 1983) and (Mr, 1992). 

It is useful to bear in mind that the real and imaginary parts of Ao,2 are always taken with respect to the CP- 
violating phase and not the final-state strong interaction phases that have already been extracted in eq. (1.9). A 
simpler form of eq. (1.11), in which ImAo = 0, is found in those papers that follow the Wu-Yang phase convention. 
In this case e = e. 

In the standard model, e' can be in principle different from zero because the 3x3 Cabibbo-Kobayashi-Maskawa 
(CKM) matrix Vij, which appear in the weak charged currents of the quark mass eigenstate, can be in general 
complex (Kobayashi and Maskawa, 1973): 
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Vud Vus Vuh\ I 1-AV2 A A\^{p-hi)\ 

Vcd Vcs VcbX^y -A - lA^\''^ 1 - AV2 . (1.16) 

Vtd Vts Vtb J \AX^{1- p-iTj) -AX^l+iX^T^) 1 / 

In eq. (1.16) we have used the Wolfenstein parameterization in terms of four parameters: A, A, rj and p and retained 
all imaginary terms for which unitarity is achieved up to O(A^), with A = \Vus\ — 0.22. On the other hand, in 
other models like the superweak theory (Wolfenstein, 1964), the only source of CP violation resides in the K'^-K^ 
oscillation, and e' vanishes. It is therefore of great importance to establish the experimental value of e' and discuss 
its theoretical predictions within the standard model and beyond. 



1. A Brief History 



The presence in nature of indirect CP violation is an experimentally well established result (Barnctt et al., 1996) 

\e\ = (2.266 ± 0.017) x 10"^ , (1.17) 

which can be understood both qualitatively and quantitatively in the framework of the standard model of electroweak 
interactions with three generations of quarks. On the other hand, after 34 years from the discovery of Christenson et 
al. there is still no conclusive experimental evidence for a non- vanishing e' . 
The ratio e' / e is measured from 

-1 + 6 Re-. (1.18) 

e 

As discussed above, a non- vanishing e'/e gives the experimental evidence for direct CP violation. Due to the accuracy 
in the counting of Kl,s decays required by the expected smallness of |e'/^l in the standard model its detection 
represents a hard experimental challenge. 
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FIG. 1. 25 years of experiments on e'/e (in units of 10 ^). The last mark on the right, at the average central value of the 
1992-93 experiments, shows the experimental precision expected in the forthcoming years. 

As it shown in Fig. 1, the experimental error in the determination of this quantity has been dramatically reduced 
over the years from lO^'^ in the 70's (Holder et al., 1972; Banner et al., 1972; Christenson et al., 1979a; Christenson 
et al., 19796) to 3.5 x 10"^ in 1985 (Black et al., 1985; Bernstein et al., 1985) and to roughly 7 x 10"^ in the last run 
of experiments in 1992 at CERN and FNAL that obtained respectively (Barr et al., 1993; Gibbons et al., 1997) 

Re e'/e = (23 ± 3.6 ± 5.4) x 10"^ (NA31) , (1.19) 
Re e'/e = (7.4 ± 5.2 ± 2.9) x 10"'' (E731) , (1.20) 

where the first error is statistical and the second one systematic. As the reader can see, the agreement between the 
two experiments is within two standard deviations. Moreover, only the CERN result is definitely different from zero. 
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Before the end of 1999 the new FNAL (E832-KTeV) (O'Dell, 1997) and CERN (NA48) (Holder, 1997) experiments 
should provide data with a precision of (1 2) x 10~"^ and hopefully settle the issue of whether e' je is or is not zero. 
Results of the same precision should also be achieved at DA$NE (KLOE) (Patera, 1997), the Prascati ^-factory. For 
a detailed account of the experimental setups and a critical discussion of the issues involved see the review article 
by (Winstein and Wolfenstein, 1993). 

Prom the theoretical point of view, the prediction of the value of e'/e has gone through almost twenty years of 
increasingly more accurate analyses. By the end of the 70's, it had been recognized that within the standard model 
with three generations of quarks, direct CP violation is natural and therefore the model itself is distinguishable from 
the superweak model. This understanding was the result of an intensive work leading to the identification of the 
dominant operators responsible for the transition, the so-called penguin operators, and the role played by QCD in 
their generation (Vainshtein et al., 1975; Vainshtein et al., 1977). Typical estimates during this period gave e'/e 
~ 10"^ - 10-2 (Ellis et al., 1976; Ellis et al., 1977; Oilman and Wise, 1979a; Oilman and Wise, 19796). 

The next step came in the 80's as the gluon penguin operators above were joined by the electromagnetic operators 
together with other isospin breaking corrections (Bijnens and Wise, 1984; Donoghue et al., 1986; Buras and Gerard, 
1987; Lusignoli, 1989). It was then recognized that these contributions tend to make e' smaller because they have 
the opposite sign compared to the gluonic penguin contributions. This part of the computation became particularly 
critical when by the end of the decade it was realized that the increasingly large mass of the t quark would lead to an 
increasingly large contribution of the electroweak penguins (Flynn and Randall, 1989; Buchalla et al., 1990; Paschos 
and Wu, 1991; Lusignoli et al., 1992). This meant a potentially vanishing value for e'/e because of the destructive 
interference between the two contributions. 

By the 90's the entire subject was mature for a systematic exploration as the short-distance part was brought under 
control by the next-to-leading order (NLO) determination of the Wilson coefficients of all relevant operators (Buras 
et al., 1992; Buras, Jamin and Lautenbacher, 1993a; Buras, Jamin, Lautenbacher and Weisz, 1993; Buras, Jamin and 
Lautenbacher, 19936; Ciuchini et al., 1993; Ciuchini et al., 1994). This theoretical achievement together with the 
discovery of the t quark (and the determination of its mass (Barnett et al., 1996)) removed two of the largest sources 
of uncertainty in the prediction. At the same time, independent efforts were brought to bear on the matrix elements 
estimate. All combined improvements made possible the current predictions of the value of e' je within the standard 
model (Heinrich et al., 1992; Paschos, 1996; Buras, Jamin and Lautenbacher, 19936; Buras et al., 1996; Ciuchini et al., 
1993; Ciuchini et al., 1995; Ciuchini, 1997; Bertolini et al., 1996; Bertolini et al., 19986) that we are to going to review. 

2. Outline 

The analysis of e'/e can be divided into the short-distance (perturbative) part and the long-distance (mainly non- 
perturbative) part. As already mentioned, the short-distance part is by now known at the NLO level and is therefore 
under control. This part of the computation is briefly reviewed in the next section. The long-distance component 
has been studied by a variety of approaches — lattice QCD, phenomenological estimates and QCD-like models — all of 
which are eventually combined with chiral perturbation theory. As the long-distance part is the most uncertain, we 
will spend most of the review on that issue. Section II and III set the common ground on which all approaches are 
based. Section IV reviews the various determinations of the hadronic matrix elements. After a brief detour, in section 
V, to determine the relevant CKM matrix elements, in sections VI and VII we bring all elements together to discuss 
some simple models. Wc then summarize the current theoretical predictions in the standard model and comment on 
the issue of new physics. 

For a broader view on CP violation which complements the present review, especially in the attention to the 
experimental issues, the reader is encouraged to consult the article previously published in this journal (Winstein and 
Wolfenstein, 1993). 

II. THE QUARK EFFECTIVE LAGRANGIAN AND THE NLO WILSON COEFFICIENTS 

The study of kaon decays within the standard model is made complicated by the huge scale differences involved. 
Energies as far apart as the mass of the t quark and the mass of the pion must be included. The most satisfactory 
framework for dealing with physical systems defined across different energy scales is that of effective theories (Weinberg, 
1980; Georgi, 1984). The transition amplitudes of an effective theory are assumed to be factorizable in high- and 
low-energy parts. The degrees of freedom at the higher scales are step-by-step integrated out, retaining only the 
effective operators made of the lighter degrees of freedom. The short-distance physics, obtained from integrating out 
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the heavy scales, is encoded in the Wilson coefficients that multiply the effective operators. Their evolution with the 
energy scale is described by the renormalization group equations (Wilson, 1971). 

Figure 2 shows the typical diagrams that in the standard model generate the operators of the effective AS' = 1 
Lagrangian. 




(a) 



(b) 



(c) 




(d) 



w 

(e) 



FIG. 2. Standard model contributions to the matching of the quark operators in the effective flavor-changing Lagrangian. 



The AS* = 1 quark effective lagrangian at a scale /x < rric can be written (Shifman et al., 1977; Oilman and Wise, 
1979a; Bijnens and Wise, 1984; Lusignoh, 1989) as 

(2.1) 



'^w = Ciifi) Qi{n) , 



where 



G 



(2.2) 



Here if the Fermi coupling, the functions Zi{ij) and yi(/i) are the Wilson coefficients and Vij the CKM matrix 
elements; t ~ ^ytdYts/^udVus- According to the standard parameterization of the CKM matrix, in order to determine 
e'/e, we only need to consider the 2/i(/x) components, which control the CP- violating part of the Lagrangian. The 
coefficients yiifj), and Zi{iJb) contains all the dependence of short-distance physics, and depend on the t, W, b, c masses, 
the intrinsic QCD scale Aqcd and the renormalization scale /i. 

The Qi are the effective four-quark operators obtained in the standard model by integrating out the vector bosons 
and the heavy quarks t, b and c. A convenient and by now standard basis includes the following ten operators: 



Qi 


= {SaUp)y 


_A ("/3C?a)v-A ' 


Q2 


= (su)v_A 


(mc?)v-a ' 




= (Srf)v-A 




Q4,6 


= {Sadl3)y. 


-A T,qiQ0Qa)vTA , 




= l{sd)y_ 


A Eg (99)v±A ' 


Q8,W 


= i(Sarf/3) 


V-A Eg 69(9/390 )v±A 



(2.3) 



where a, /3 denote color indices {a,(3 = 1, . . . , Nc) and are the quark charges (e„ = 2/3, ed = is = —1/3). Color 
indices for the color singlet operators are omitted. The labels {V ± A) refer to the Dirac structure 7^(1 ± 75). 

The various operators originate from different diagrams of the fundamental theory. First, at the tree level, we 
only have the current-current operator Q2 induced by Vl^-exchange. Switching on QCD, a one-loop correction to 
Ty-exchange (like in Fig. 2b, c) will induce Qi. Furthermore, QCD through the penguin loop (Fig. 2d) induces the 
gluon penguin operators Qs-e- The gluon penguin contribution is split in four components because of the splitting 
of the gluonic coupling into a right- and a left-handed part and the use of the SU {Nc) relation 
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2 T^gT^is = Sa/iS^s - j^SasSj/s , (2.4) 

where is the number of colors, a = 1, — 1 and T° are the properly normalized SU{Nc) generators, 

TrTT*" = 1/2 (5"'', in the fundamental representation. Electroweak loop diagrams — where the penguin gluon is 
replaced by a photon or a Z-boson and also box- like diagrams — induce Qr^g and also a part of Qs. The operators 
Qs,io are induced by the QCD renormalization of the electroweak loop operators Qr^g. 

Even though the operators in eq. (2.3) are not all independent, this basis is of particular interest for any numerical 
analysis because it is that employed for the calculation of the Wilson coefficients to the NLO order in as and ae (Buras 
et al., 1992; Buras, Jamin and Lautenbacher, 1993a; Buras, Jamin, Lautenbacher and Weisz, 1993; Buras, Jamin and 
Lautenbacher, 19936; Ciuchini et al., 1993; Ciuchini et al., 1994) and we will use it throughout. 

Anticipating our discussion, the pie chart in Fig. 3 shows pictorially the relative importance of the operators in 
eq. (2.3) in the final determination of the value of e'/e , as obtained in the vacuum saturation approximation to 
the hadronic matrix elements. In particular, Fig. 3 shows the crucial competition between gluonic and electroweak 
penguins in the determination of the value of e'/e . Such a destructive interference might accidentally lead to a 
vanishing e'/e even in the presence of a source of direct CP violation. This feature adds to the theoretical challenge 
of predicting s' /e with the required accuracy. 




FIG. 3. Relative contributions to e'/e of the operators in eq. (2.3). Operators giving a (negative) positive contribution are 
depicted in (dark) light gray. All matrix elements are taken in the vacuum saturation approximation. The operators Qs and 
Qio are generated by QCD and do not receive contributions from the one-loop matching. 

While there exist other possible operators in addition to those listed in eq. (2.3), they are numerically irrelevant 
within the standard model. For instance, the two operators 

Qii = -At s [uidR + msL] a ■ G d and Q12 = — ^ s [rudR -f rUsL] a ■ F d , (2.5) 

where i? = (1 + 75) /2 and L = {1 ~ 75) /2, are present. These operators are induced by gluon and photon penguins 
with a free gluon (photon) leg. The matrix elements of these operators give a vanishingly small contribution to 
K ^ TTTT decays (Bertolini et al., 1995; Bertolini et al., 1998a). 

In table I we summarize in a synthetic way the diagrammatic origin of the contributions to the various Wilson 
coefficients when considering the one-loop matching of the quark effective lagrangian with the full electroweak theory. 
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TABLE I. Contributions to the one- loop matching of the AS = 1 Wilson coefficients at /^ = mw- The notation refers to that 
of cq. (2.1) and Fig. 2. Non-vanishing contributions to Cs and Cio arise via the QCD renormalization of the operators Qr and 
Qg, respectively. 



/i — rriw Ci C2 C3 C4 C5 Ce C7 Cs Cg C] 

Tree 7 

Tree + g ^ \j 
Tree 7 V 

^3 V V V 7 



Having established the operator basis, a full two-loop calculation (up to and asap-m) of the quark operator 
anomalous dimensions is performed. This calculation allows us via renormalization group methods to evaluate the 
Wilson coefficients at the typical scale of the process, thus resumming (perturbatively) potentially large logarithmic 
effects to a few 10% uncertainty. As already mentioned, the size of the Wilson coefficients at the hadronic scale (of 
the order of 1 GeV) depends on and the threshold masses mt, mw ■, rub and me- The top quark mass dependence 
enters in the penguin coefficients yi{^i) via the initial matching conditions for the renormalization group equations. 

Small differences in the short-distance input parameters are present in the various treatment in the literature. In 
order to give the reader an idea of the ranges used, we list below some of the values. 

The most recent determination of the running strong coupling in the MS scheme is (Barnett et al., 1996) 

a^(TOz) =0.119 ±0.002, (2.6) 

which at the NLO corresponds to 

A^^j5 = 340 ± 40 MeV . (2.7) 

For mt we take the value (Tipton, 1997) 

mr''' = 175±6 GeV. (2.8) 

The knowledge of the top quark mass is one important ingredient in the reduced uncertainty of the recent estimates 

of s'/s . 

The relation between the pole mass M and the MS running mass m(/i) is given at one loop in QCD by 



m(M) = M(g2 = M^) 



^ 4 q,(M) 

3 TT 



(2.9) 



For the running top quark mass, in the range of ag considered, we then obtain 

mt{mt) ~ 167 ± 6 GeV , (2.10) 

which, using the one-loop running, corresponds to 

mt{mw) ^Vn^l GeV , (2.11) 

which is the value to be used as input at the mw scale for the NLO evolution of the Wilson coefficients. In eq. (2.11) 
we have averaged over the range of Aq^j-, given in eq. (2.7). 

The use of the running top mass in the initial matching of the Wilson coefficients softens the matching scale 
dependence present in the LO analysis. By taking [i = ?7if°''' as the starting matching scale in place of mw, and using 
correspondingly mt{mt), the NLO Wilson coefficients of the electroweak and gluon penguins at /z ~ 1 GeV, remain 
stable up to the 10% percent level. 

For mt we have the mass range (Barnett et al., 1996) 
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mP°''' = 4.5 H- 4.9 GeV , 



(2.12) 



which corresponds to 

mb{mb) = 4.1 ^ 4.5 GeV . (2.13) 

Analogously, for rric one has 

mP°''' = 1.2 H- 1.9 GeV, (2.14) 

which corresponds to 

mc(mc) = 1.0 ^ 1.6 GeV , (2.15) 

Values within the MS ranges have to be used as the quark thresholds in evolving the Wilson coefficients down to the 
low-energy scale where the matching with the hadronic matrix elements is to be performed. 

In choosing the quark mass thresholds one should bear in mind that varying ml°^'^ within the given range affects the 
final values of the Wilson coefficients only at the percent level, while varying the charm pole mass in the whole range 
given may affects the real part of the gluon penguin coefficients up to the 20% level. We will take mb{mb) = GeV 
and mc(mc) = 1.4 GeV. 

TABLE II. The AS = 1 NLO Wilson coefficients relevant for CP violation arc given at /i = 1.0 GeV for mt{mw) = 177 GeV, 
which corresponds to m^"^" = 175 GeV, (a = 1/128). In addition one has yi,2{n) = 0. 



^(4) 


300 MeV 






340 MeV 


380 MeV 






HV 


NDR 


HV 


NDR 


HV 


NDR 


ys 


0.0341 


0.0298 


0.0378 


0.0326 


0.0420 


0.0356 


2/4 


-0.0558 


-0.0530 


-0.0597 


-0.0564 


-0.0639 


-0.0597 


2/s 


0.0149 


0.000687 


0.0160 


-0.00204 


0.0173 


-0.00581 


2/6 


-0.0883 


-0.100 


-0.0994 


-0.115 


-0.113 


-0.133 


yr/a 


-0.0202 


-0.0210 


-0.0195 


-0.0209 


-0.0188 


-0.0209 


ys/a 


0.184 


0.169 


0.209 


0.192 


0.240 


0.220 


yg/a 


-1.70 


-1.70 


-1.75 


-1.74 


-1.80 


-1.80 


2/10 /a 


0.735 


0.722 


0.806 


0.790 


0.885 


0.867 



In table II we report the numerical values of the NLO Wilson coefficients relevant for CP violation in A5 = 1 

processes. The coefficients yi(/i) are given at the scale /i = 1 GeV and are dependent on the choice of the 75 
scheme in dimensional regularization. The values in the table refer to two commonly used schemes, namely the naive 
dimensional regularization (NDR), in which 75 anticommutes with the Dirac matrices in d dimensions, and the t' 
Hooft-Veltman scheme (HV) ('t Hooft and Veltman, 1972), in which they anticommute only in four dimensions. The 
latter prescription has been shown to be a consistent formulation of dimensional regularization in the presence of 
chiral couplings (Breitenlohner and Maison, 1977). A consistent calculation of the hadronic matrix elements should 
match the unphysical scale and scheme dependence of the Wilson coefficients so as to produce a stable amplitude at 
the given order in perturbation theory. We will return on this issue in Sect IV when discussing the various approaches 
to the long distance part of the calculation. 

The case of the AS = 2 theory is treated along similar lines. The effective AS = 2 quark lagrangian at scales 
/U < nic is given by 

£as=2 = -C2s(m) Qs2{fi) , (2.16) 

where 

C2si^^) = [>^lviS{xc) + \lmS(,xt) + 2X,Xtrj3Sixe, Xt)] h{y) (2.17) 

where Xj = Vj^V^* , Xi = mf/m^. We denote by Qs2 the AS = 2 local four quark operator 
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QS2 = {SL7^dL){SL7ij.dL) , 



(2.18) 



which is the only local operator of dimension six in the standard model. 

The integration of the electroweak loops leads to the Inami-Lim functions (Inami and Lim, 1981) S{x) and S{xc, Xt), 
the exact expressions of which can be found in the reference quoted, depend on the masses of the charm and top 
quarks and describe the AS* = 2 transition amplitude in the absence of strong interactions. 

The short-distance QCD corrections are encoded in the coefficients rji, r]2 and 773 with a common scale-dependent 
factor b{ii) factorized out. They are functions of the heavy quarks masses and of the scale parameter Aqcd- These 
QCD corrections are available at the NLO (Buras et al., 1990; Herrlich and Nierste, 1994; Herrlich and Nierste, 
1995; Herrlich and Nierste, 1996) in the strong and electromagnetic couplings. 

The scale-dependent common factor of the short-distance corrections is given by 

2/9 f. J (m) 



6(m) = [as (m)]-'^' [1 - ) , (2.19) 

where J3 depends on the 75-scheme used in the regularization. The NDR and HV scheme yield, respectively: 



,Noa_307 .r = -^. (2.20) 



All the other numerical inputs can be taken as in the AS' = 1 case. 



III. CHIRAL PERTURBATION THEORY 



Quarks are the fundamental hadronic matter. However, the particles we observe arc those built out of them: 
baryons and mesons. In the sector of the lowest mass pseudoscalar mesons (the would-be Goldstonc bosons: tt, K and 
jy), the interactions can be described in terms of an effective theory, the chiral lagrangian, that includes only these 
states. The chiral lagrangian and chiral perturbation theory (Weinberg, 1979; Gasser and Leutwyler, 1985; Gasser 
and Leutwyler, 1984) provide a faithful representation of this sector of the standard model after the quark and gluon 
degrees of freedom have been integrated out. The form of this effective field theory and all its possible terms are 
determined by SULi'i) x SUii{3) chiral invariance and Lorentz invariance. The parts of the lagrangian which explicitly 
break chiral invariance are introduced in terms of the quark mass matrix M. 

The strong chiral lagrangian is completely fixed to the leading order in momenta by symmetry requirements and 
the Goldstone boson's decay constant /: 



,(2) 



■^strong ^ 

where M. = diag[m«, nid, rus] and Bq is given by {(jiqj) 



Tr (D^T.D'^T,^) + ^—Bq Tr {MY:^ + T,M^) 



2 

pBoSij, with 



(3.1) 



K 



mu + md 



(3.2) 



according to PCAC in the hmit of SUi'S) flavor symmetry (fK = U = /)• The SUl{S) x SUr{3) field 



'2i 

exp ( -J U{x) 



(3.3) 



contains the pseudoscalar octet: 



1 ^ 1 
n(x) = - V A»7r"(x) = 

a=l * 



Tr° 7r+ K+ 
K- ^0 ^8 



(3.4) 



where 



7T° = 



V2 \/6 



V2 \/6 



\/6 



(3.5) 
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The coupling / is, to lowest order, identified with the pion decay constant /,r (and equal to //f before chiral loops are 
introduced); it defines a characteristic scale 



= 27rv/677Vc/ ~ O.SGeV, (3.6) 

typical of the vector meson masses induced by the spontaneous breaking of chiral symmetry. When the matrix S is 
expanded in powers of /~^, the zeroth order term is the free Klein-Gordon lagrangian for the pseudoscalar particles. 
Under the action of the elements Vr and Vl of the chiral group SUii{3) x SUl{3), S transforms linearly: 

E' = Vr1:vI , (3.7) 

with the quark fields transforming as 

q'l = Vl Ql and q'^ = VRqR , (3.8) 

and accordingly for the conjugated fields. 

Quark operators are represented in this language in terms of the effective field S and its derivatives. For instance, 
at the leading order, the quark currents are given by 



-'--i^i^^D.E)^^, (3.9) 



2 

f2 



while the quark densities can be written at O(p^) as 



'^~"'l^-^-^(ED^St).. , (3.10) 



^ -2Bo 



f2 



4 



^Et + is Et£)^E£>''E+ + ABqLsT^^M E^ 



(3.11) 



where Lz,^ are coefficients which belong to the 0{p'^) chiral lagrangian. To the next-to-leading order in the momenta, 
in addition to the leading order chiral lagrangian (3.1), there are ten chiral terms and thereby ten coefficients Li to 
be determined (Gasser and Lcutwylcr, 1984; Gasscr and Lentwyler, 1985) either experimentally or by means of some 
model. As we shall see, some of them play an important role in the physics of e'/e . As an example, we display the 
is and Z/8 terms in '^strong which, appecir in ec[. (3.11) cind govern much of the penguin physics: 

is Bo Tr [il^Eil^Et (XE^ + EA4^)] and ig Bo Tr [M^EM^T, + M'S^MT,^ . (3.12) 



A. The Weak Chiral Lagrangian 

We can write the most general expression for the AS' = 1 chiral lagrangian in accordance with the SU{3) l x SU{3)r 
symmetry, involving unknown constants of order Gf- This is done order by order in the chiral expansion. Typical 
terms to 0{p^) are obtained by inserting appropriate combinations of Gell-Mann matrices into the strong lagrangian. 
The corresponding chiral coefficients must then be determined by means of some model or by comparison to the 
experimental data. 

We find it convenient to write the AS = 1 chiral lagrangian at 0{p^) in terms of the following eight terms, of which 
seven are linearly independent: 

4'Li= Gi°],(Qr,8)Tr(A3EtA}E) 

+ G^LRiQ7.s) [Tr(A^EU}SM^S) +Tr(AlEAiEl'A^Et)] 

+ Gs{Q3-io)Tt{XID^J:^D^^) 

+ G2i(Qi,2,9,io) Tr (A?Eti?,,E) Tr (A^Et^^E) 

+ Gi^(Qi,2,9,io) Tr (A^Eti^^^E) Tr {XIj:W'^J:) 

+ GlRiQ7,8) Tr (A^T^^EA^i^^Et) 

+ GiRiQrs) Tr (A^Et^^E) Tr (A}Ei?"Et) 

+ GIr{Q7,8) [Tr (A?E) Tr (A^i)^Eti>''E Et) + Tr (A?i>^Ei)''Et E) Tr (A^Et)] , (3.13) 
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where A* are combinations of GcU-Mann SU{3) matrices defined by {/(^j)ik = SuSjk and S is defined in cq. (3.3). The 
covariant derivatives in eq. (3.13) are taken with respect to the external gauge fields whenever they are present. Other 
terms are possible, but they can be reduced to these by means of trace identities. 

The non-standard form and notation of cq. (3.13) is chosen to remind us of the flavor and chiral structure of the 
effective four-quark operators which are represented by the various terms. In particular, in Gs we collect the (8^, x Ij^) 
part of the interaction which is induced by the gluonic penguins and by the analogous components of the electroweak 
operators Qj-io- The two terms proportional to G^/^ and G^^j^ arc an admixture of the (27^ x 1^) and the (8^^ x ]_^) 
part of the interactions induced by the left-handed current-current operators Qi, 2,9.10- The term proportional to 
is the constant (non-derivative) 0(p") part arising from the isospin violating (8^ x 8^) electroweak operators. The 
0(p^) corrections to Gj°^ are the quark mass term proportional to G^"^^ (related to L^), the momentum corrections 
proportional to Glj^ (related to L5) and G^']^. One may verify that G^-^ and G^j^ can be obtained by multiplying the 
bosonized expression of a left- and a right-handed quark density (in a manner similar to Qq), while G^^ is obtained as 
the product of a left- and a right-handed quark current. It is therefore natural to call these terms factorizable (although 
GJ^^ has a non-factorizable contribution). The term GJ^ is, however, genuinely non-factorizable (Fabbrichesi and 
Lashin, 1996). 

The terms proportional to Gs, GJ^ and G^j^ have been studied in the literature (Cronin, 1967; Pich and de Rafael, 

1991; Bijncns ct al., 1993; Eckcr ct al.. 1993) in the framework of chiral perturbation theory. The three terms are not 
independent. Those proportional to G^^ and G'l^^ can be written in terms of the 8 and 27 SU{3)l components as 
follows: 



£27 = G27{Qi) 



^ Tr (A?Et£)M5]^ xr (A^E^'D^S) ^ Tr (A^E+D^E) Tr (AJE^D^E) 



(3.14) 



3 

which transforms as (27^;, x 1^), and 

C-i = G^iQC) [Tr (AjEt/^'^E) Tr (A^St^^E) - Tr (AiEt£>^E) Tr (AjE^i^^E)] , (3.15) 

which transforms as (8^ x 1^). We prefer to keep the AS* = 1 chiral Lagrangian in the form given in eq. (3.13), which 
makes the bosonization of each quark operator more transparent, and perform the needed isospin projections at the 
level of the matrix elements. Equations (3.14)-(3.15) provide anyhow the comparison to the standard notation. The 
chiral coefficients in the two bases are related by 

G8(Qi) = i [3 G2 JQO - 2 Gli(QO] 

Grim = I [Gli^iQi) + GL(Qi)] , (3.16) 

for i = 1,2. Notice that there is no over-counting of the 8^ x l^j contributions to eq. (3.13) from the operators Qg.io 
when a consistent prescription like that given in (Antonelli et al., 1996) is followed. 

Concerning the (8^ x 8^) part of the AS* = 1 chiral lagrangian, the constant term was first considered in (Bijncns and 
Wise, 1984), while its mass and 0{p^) momentum corrections were first discussed in (Antonelli et al., 1996; Bertolini 
et al., 19986). 

As an example of the form of the chiral coefficients, we give the determination in the leading order in 1 /N^. of the 
two most important contributions to e'/e : 

G8(Q6) = -24M!^C'6 (3.17) 

and 

G^'\Qs) = -3{qqf Gs , (3.18) 

where Ge.s are the Wilson coefficients of the operators Qq.s at the matching scale /i. 

The AS' = 1 O(p^) Lagrangian is much more complicated (Kambor et al., 1990; Esposito-Farese, 1991; Ecker et al., 
1993; Bijnens et al., 1998) but we will not need its explicit form. In fact, only certain combinations of coefficients 
from the O(p^) arc required in order to compute the relevant amplitudes to this approximation. 

The AS = 2 weak chiral lagrangian is simpler. At the leading order 0{p^), the AS = 2 weak chiral lagrangian is 
given by only one term: 
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4's=2 = G{Qs2) TV (A^ED^Et) TV (A^SD'^St) . (3.19) 
The chiral coefficient is in this case given at the LO in l/Nc by 

G{Qs2) = -^-^C2s- (3.20) 



IV. HADRONIC MATRIX ELEMENTS 



The estimate of the hadronic matrix elements must rely on long-distance effects of QCD. It is useful to encode the 
result of different estimates in terms of the -B, parameters that are defined in terms of the matrix elements 

(Oi)o,2 = ((7r^)(/=o,2)|Q.|i^°) (4.1) 



as 



p(o,2) _ J^^w»;o,2 

= \VSA ' 



{Q^ 



0,2 



and give the ratios between hadronic matrix elements in a model and those of the vacuum saturation approximation 
(VSA). The latter is defined by factorizing the four-quark operators, inserting the vacuum state in all possible manners 
(Ficrzing of the operators included) and by keeping the first non- vanishing term in the momentum expansion of each 
contribution. 

As a typical example, the matrix element of Qe in the factorized version can be written as the product of density 
matrix elements 

(7r+7r-|Q6|^°) = 2(7r+|u75d|0)(7r-|su|if") -2(7r+7r-|dd|0)(0|s75rf|/sr°) 

+2 [(0|ss|0) - {0\dd\0)] (7r+7r-|s75rf|iV'°), (4.3) 

where the matrix elements like (0|s75u|if+) and (tt+I |_ftr+) are obtained from PCAC and the standard parame- 
terization of the corresponding currents, (0| S7'' (1 — 75) u \K~^) and (tt+I S7'' (1 — 75) u \K~^). In the same way, the 
left-left currents operators can be written in the factorizable approximation in terms of matrix elements of the currents. 

Notice that the definition in eq. (4.2) neglects the imaginary (absorptive) parts of the hadronic matrix elements. 
Imaginary and real components, when multiplied by the corresponding short-distance coefficients and summed over 
the contributing operators, should reproduce the global phase of the amplitude arising from final state interactions. 
However, some approaches to hadronic matrix elements do not account for absorptive contributions. Therefore, 
in order to make the discussion of the Bj factors of different models as homogeneous as possible, we propose the 
definition in eq. (4.2). Consistently with the use of such a definition, extra overall l/cos(5o,2 factors appear in the 
7 = 0, 2 amplitudes, as discussed in Sect. VI. 



A. Preliminary Remarks 

The Bi parameters depend in principle on the renormalization scale /x and therefore they should be given together 
with the scale at which they are evaluated. 

In this respect, in a truly consistent calculation of the hadronic matrix elements, the cancellation of the unphysical 
renormalization scale and scheme dependence of the Wilson coefficients should formally be proven order by order in 

perturbation theory. 

The only approach that fully satisfies these requirements is that based on the lattice regularization (discussed in 
subsection F), where the same theory, namely QCD, is used in both the short- and the long-distance regimes and the 
matching only involves the different regularization schemes. 

The Miinchen phenomenological approach (discussed in subsection E) represents a clever attempt to address the 
problem of a consistent calculation of s'/s in a framework originally based on the 1/A^c expansion. In this approach 
one extracts as much information as possible on the hadronic matrix elements by fitting the AJ =1/2 selection rule 
at a fixed scale and in a given renormalization scheme. The scale and renormalization scheme stability of physical 
amplitudes can then be obtained using perturbation theory since the matching scale between short- and long- distance 
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calculations is large enough (/i = rric) to lie inside the perturbative regime. The phenomenological input allows for a 
direct determination of the current-current matrix elements and indirectly of some of the penguin matrix elements, 
thus reducing the number of free parameters in the AS = 1 effective lagrangian. On the other hand, the same fit 
does not give any information on the actual value (and scheme dependence) of the Be, 8 parameters at the given scale, 
which are the most relevant for determining e'/e . 

In the Trieste group approach (discussed in subsection G) there is no attempt to prove formally the consistency 
of the matching along the lines stated above. The matching is done between QCD on the short-distance side and 
phenomenological models, the xQM and chiral perturbation theory, on the long-distance side. In the long-distance 
calculation the scale and renormalization scheme dependences appear naturally. It is then assumed that these un- 
physical dependences may satisfactorily match those of the short-distance calculation. The fact that this assumption 
is numerically verified (even beyond expectation), thus giving at the given order of the calculation a stable set of pre- 
dictions, and that it allows for a complete calculation of all matrix elements in terms of a few basic "non-perturbative" 
parameters, make this phenomenological analysis valuable. The pattern of contributions which emerges and which 
leads to a satisfactory reproduction of the AI =1/2 rule may be of help in other investigations. The major weakness 
of the approach is the poor convergence of the chiral expansion at matching scales of the order of the p mass or higher, 
which are required by the reliability of the perturbative strong coupling expansion. 

Very recently the Dortmund group (see subsection D) has developed a systematic procedure for matching short- and 
long-distance calculations, improving both technically and conceptually on the original 1/Nc approach of (Bardeen 
et al., 1987). On the other hand, at the present status of the calculation, the scale stability of the matching with 
the short-distance coefficients is for some of the relevant observables (A/ = 1/2, 3/2 amplitudes, Bk) quite poor 
(Hambye, 1997; Kohler, 1998). 

B. The Vacuum Saturation Approximation 

According to the discussion above it is clear that there is no theoretical underpinning for the consistency of the 
VSA; it is a convenient reference frame which is equivalent to retaining terms of 0(1/ N^^ in the l/A'^c-expansion to 
the leading (non-vanishing) order in the momenta for all Fierzed forms of the operators. Its application should in 
general not be pushed beyond leading order in the strong coupling expansion. On the other hand, we find it useful 
for illustrative purposes to use the VSA hadronic matrix elements together with NLO Wilson coefficients in order 
to exhibit some features of the long-distance calculation and allow for a homogeneous comparison with the other 
estimates. For this purpose we will use in all numerical estimates the Wilson coefficients obtained in the HV scheme 
and set the matching scale at 1 GcV (see table II). 

Some of the relevant VSA hadronic matrix elements depend on parameters that are not precisely known. As a 
consequence, the knowledge of the Bi is not the whole story and, depending on assumptions, different predictions 
of e'/e may well differ even starting from the same set of Bi. It is therefore important to define carefully the VSA 
matrix elements. According to the standard bosonization of currents and densities at O(p^) one obtains: 



{Ql)2 



{Q2)2 



{Qi)o 



(Q3)0 



(Q2)o 



(Q4)o 




16 {qq)^Lr, 



(4.6) 



(4.7) 



(4.4) 



(4.5) 



(4.8) 
(4.9) 



(Q5)0 



X 



(4.10) 



((96)0 




(4.11) 
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(Q7)0 
(07)2 

{Qs)o 

{Qs)2 
(Q9)0 
{Q9)2 

(Qio)o 

(Ql0)2 



2^/3 {qqf , 8 



2 



X 



-X 



p 



X 



2Nr. 



-X 



V2 



-X 



2' 

-I- 



1 + 



2N, 
1 

1 



1 



1 



where 



X = V3f{m%-ml) . 

In addition, from the O(p^) chiral lagrangian evaluation of fxj f-K one obtains, neglecting chiral loops, 



1 (J k-P 



p 



m 



K 



mi 



(4.12) 
(4.13) 
(4.14) 
(4.15) 
(4.16) 
(4.17) 
(4.18) 
(4.19) 

(4.20) 

(4.21) 



while the quark condensate may be written in terms of the meson and quark masses using eq. (3.2). The subleading 

1/Nc terms arise from the Fierzing of the quark operators via the SU{Nc) relation (2.4). 

In a similar manner, in the case of the AS = 2 amplitude, the scale-dependent Bk parameter is defined by the 
matrix element 



The scale independent parameter Bk is defined by 

BK = b{^x)BK{^l). 

In the VSA, for which 6(/x) = 1, the value 



Bk = 



1 



1 



(4.22) 



(4.23) 



(4.24) 



is found. 

As it has been mentioned before, already at the level of the VSA, it is necessary to know the value of /, (qq) 
or, via PCAC, the value of quark masses. Specifically, unless otherwise stated, we will assume as reference values 
for the input parameters in the VSA / = and {qq) (1 GeV) = —(238 MeV)"^, which corresponds via eq. (3.2) to 
(m„ + md){l GeV) = 12 MeV, or equivalently to (m^ + md)(l GeV) = 157 MeV. 

Notice that the evaluation of the matrix elements of the operators Qe-s requires already at the VSA level the 
strong 0{p'^) chiral coefficient L5. For this reason, the determination of Bq has been disputed in the past (Dupont 
and Pham, 1984; Gavela et al., 1984; Donoghue, 1984; Chivukula et al., 1986). 

We shall discuss the numerical results of the Bi factors in an improved VSA model which includes the complete 
0{p^) corrections to the leading momentum independent terms in the Qr,8 matrix elements. In the same model we 
will show the effect of the inclusion of final state interactions. Then, we will summarize the published results of the 
three most developed estimates: the Miinchen phenomenological approach, the Roma numerical simulations on the 
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lattice and, among possible effective quark models, the chiral quark model (for which the complete set of operator 
basis has been analyzed by the Trieste group). 

The values quoted for the Bi are taken at different scales so that they cannot be directly compared. Notice, however, 

(2) 

the two most important parameters, namely Bq and Bg ' have been shown to depend weakly on the renormalization 
scale for ;u ^ 1 GeV (Buras, Jamin and Lautenbacher, 19936). 



C. A Toy Model: VSA+ 

A comparison between the VSA matrix elements and the chiral lagrangian of eq. (3.13) shows that none of the 
0{p^) terms proportional to G'^"'-', G^j^ and G"£^ is included in the standard VSA. These contributions enter as 
additional corrections to the 0(p°) leading term in the matrix elements of the operators Qr and Qs (Antonelli et al., 
1996; Bertolini et al., 19986). With the help of eq. (3.11) and keeping all terms one obtains 

2^/3 {qqf 8 {qq)^L, ^ ^ 1 ^ ^ 16a/3 {qq)^ loT t ^ 2 
in\ -"^ V2 8V6 {qqr . . 2 u OfiA 

(Qs)o = 2 73 M! + sMj^ ^ + + 16 (2i« -L,)ml, (4.27) 

(^8)2 = ^ - + (2^8 - < ' (4-28) 

where we have neglected rn^/m^ terms. The 0(p^) wave- function renormalization has been included by multiplying 
the 0{p°) term by 

V^Z. = l-4£5 . (4.29) 

In this toy model, which we call VSA-I-, we neglect all chiral loop corrections, even though they are of O(p^) on the 
constant term in the AS* = 1 chiral lagrangian (all other chiral loop corrections are of 0(p^)). The parameter / in 
the 0(p°) terms of eqs. (4.25)-(4.28) may be rewritten in terms of the renormalized Jk and/or /,r. At O(p^) such 
a rewriting is not unique. For the purpose of the present discussion we take, as in the standard VSA, / = The 
terms proportional to 2Ls — represent the additional corrections to the VSA matrix elements. 

In order to obtain an estimate of the combination 2Lf^ — L5 consistent with that of in eq. (4.21), used in the 
VSA, we employ the mass relation (Gasser and Leutwyler, 1985) 

4 = ^^(1 + Am), (4.30) 
where m = (m„ + ma) 12 and, neglecting chiral loops, 

Am = {m\ - ml) \2U - U\ ■ (4.31) 



Assuming PCAC to hold with degenerate quark condensates, and keeping 7^ /tt, we then obtain 



2^8 ~ = — 



- 1 



f 



2 



(4.32) 



'^K '""IT 

The purpose of introducing the VSA-|- model is to show the relevance of the O(p^) corrections to the leading term 
for the (Q8)2 matrix element which is crucial in determining e'/e. The coefficients By and Bg are modified from 
their VSA values as shown in Table III. Their values are essentially independent on the value of {qq) , because of the 
smallness of the terms not proportional to the quark condensate. 
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TABLE III. The Bi in the VSA+ model described in the text. All other Bi parameters are equal to unity. 





— ^8 


0.7 




_5(2) 


0.6 



Much uncertainty in the present toy model is hidden in the approximations made in giving L5 and Lg. As an 
example, a determination of these coefficients in chiral perturbation theory including dimensionally regularized chiral 
loops gives, at the scale nip, a ' greater than one (Fabbrichesi and Lashin, 1996). 

A discussion of the implications of the VSA+ model for e' /e and a pedagogical comparison with the standard VSA 
are presented in Sect. VI. 



D. 1/Nc Corrections 

Chiral-loop corrections arc of order l/Nc and of order 0{p^) in the momenta (except for those of the leading 
electroweak term that are of Oip^)). They have been included in the 1/A^c approach of (Bardeen et al., 1987) by 
means of a cut-off regularization that is then matched to the short-distance renormalization scale between 0.6 and 1 
GcV. The values thus found = 5.2, B^^ = 2.2, ^ = 0.55) although encouraging toward an explanation of the 

A/ =1/2 rule were still unsatisfactory in view of trusting the approach for a reliable prediction of e'/e . 

Along similar lines, the Dortmund group (Heinrich et al., 1992) included chiral corrections to the relevant operators 
Qq and Qs- They did not report explicit values for their B^. However, from their analysis it is clear that they find a 
rather large enhancement of Bq and a suppression of B^. More recently (Hambye et al., 1998) have estimated these 
coefficients in a new study which pays special attention to the matching between the renormalization scale dependence 
of chiral loops, regularized by a cut-off, and the dimensionally regularized Wilson coefficients. They find almost no 
enhancement in the iJg but a larger suppression of Bg. No new calculation of e'/e has appeared so far. Some of the 
relevant observables, as Bk and the / = 0, 2 amplitudes, show at the present status of the calculation a quite poor 
scale stability (Hambye, 1997; Kohler, 1998), which may frustrate any attempt to produce a reliable estimate of e' /e . 

The parameter Bk has been independently estimated in the l/A^c expansion with explicit cut-off by (Bijnens and 
Prades, 1995), finding values between 0.6 and 0.8. 

A systematic study of chiral-loop corrections in dimensional regularization was performed first by (Kambor et al., 
1991) and more recently redone using the MS scheme by the Trieste group (Bertolini et al., 1996; Bertolini et al., 
19986). The chiral-loop corrections also generate an absorptive part in the amplitudes which should account for the 
final state interactions. In any case, they seem to play an important role in the determination of the hadronic matrix 
elements. 



E. Phenomenological Approach 

The phenomenological approach of the Miinchen group (Buras, Jamin and Lautenbacher, 1993&; Buras et al., 1996) 

writes all hadronic matrix elements in terms of just a handful of Bi: B^^ for the {V — A) (V — A) operators and 

(2) 

Be and B^ ' for {V — A) ® {V + A) operators. This approach exploits in a clever manner the available experimental 
data on the amplitudes Aq and A2 in order to extract the (scheme dependent) values of B^°2^^ and, via operatorial 

relations, of some of the penguin matrix elements, while leaving Bq and B\ ' as free input parameters to be varied 

within given limits. 

In particular, B\ 2 are obtained directly from the experimental value 

ReA2 = 1.50 X 10"^ GeV , (4.33) 
via the matching condition at fj, = rric and the scale independence of the physical amplitude as 

(gi>2 = (Q2>2 = , (4.34) 

c z+ [rric ) 
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where c = G^V^udKIs/^ ^ii^ z+ is the real part of the Wilson coefficient of the operator Qi + Q2; -Bg^io sue then 
obtained by using the operatorial relation 

(Q9,10)2 = l{Qi)2 . (4.35) 

B^l g are similarly expressed as functions of B^^ by means of other operatorial relations and matching conditions 
at the charm-mass scale. In fact, in the HV scheme at nic there are no penguin contributions to CP conserving 
amplitudes and in the NDR the penguin contamination is numerically small. Therefore one can write 

{Qi)o = 7 — r — ■z{Q2)o ■ (4.36) 

czi(mc) zi(TOc) 

Finally, b!^''' is also obtained under the plausible assumption {Q2 — Qi) > {Q2 + Qi) > 0, valid in all known non- 
perturbative approaches, from the experimental value of 

Re Ao = 33.3 x 10"* GeV . (4.37) 

The following operatorial relations, which hold exactly in the HV scheme, may then be used 

((34)0 = (Q3)o + {Q2)o - (Qi)o , (4.38) 

(<39)o = ^(Qi)o-^(<33)o, (4.39) 

(Qio)o = (Q2)o + l{Qi)o - liQsh ■ (4.40) 

It is important to recall that B^ is taken equal to 1, which may be a rather crucial assumption in the determination 
of -B4, as we shall see. 

After imposing that B5 = Bq and B^^^ = B^^ , this leaves us with only two free input parameters Bq and B^^ that 
are varied within 20% from unity. 

The parameter B^ is pragmatically taken to span from the central value of the lattice (see the next section) to 
that of QCD sum rules (Narison, 1995). 

TABLE IV. The B, in the Miinchen phenomenological approach. The results for -81,2,9,10 are obtained by fitting the A7 = 1/2 
selection rule in K ^ nn decays at the matching scale /i — rric. We show the values obtained in the HV scheme for the central 
value of Aqqq = 325 MeV. The value for B4, is obtained by assuming Bz = 1. All the remnant Bi are taken equal to 1 except 

for Bq and Sg that are varied within ±20% from unity. The parameter Bk is scale and renormalization scheme independent. 

Sf^ 13 

Bf 6.2 

> = Bf^ 0.47 

Bi 5.2 



B'^" 7.1 
SS' 7.7 
^ = B^o' 0.47 



F. Lattice Approach 

The regularization of QCD on a lattice and its numerical simulation is the most satisfactory theoretical approach 
to the computation of the hadronic matrix elements (for a review see, for instance, (Sharpe, 1994)), and should, in 
principle, lead to the most reliable estimates. However, technical difficulties still plague this approach and only some 
operators have been precisely determined on the lattice. In addition, the use of approximations like quenching make 
it very difficult to assess the effective uncertainty of the calculation. 
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Another problem of the approach is that it is stiU not possible to directly compute the K tttt amplitude in 
Euclidean space. It is therefore necessary to rely on chiral perturbation theory in order to obtain the amplitude with 
two final pions from that with just one. In this sense even the lattice approach is not, at least for the time being, a 
first-principle procedure. As a matter of fact, when considering the complete O(p^) chiral lagrangian of eq. (3.13) a 
problem arises in so far as the term proportional to G^^^ has a vanishing contribution io K ^ n. 

TABLE V. The coefficients obtained in the Roma lattice calculation at the matching scale n = 2 GeV in the NDR scheme. 
The values of B^^^ are derived from the phenomenological fit of the A7 =1/2 rule. Accordingly, Ba is varied in the range 
1-^6. All others Bi are taken equal to 1. 



1.0 ±0.2 
0.6 ±0.1 

0.8 ±0.15 
0.62 ± 0.10 
0.75 = 0.15 



Bsfi 

?(2) 
7 

(2) 



B. 



(2) 



Bk 



Table V summarizes the values obtained by direct lattice computations and used by the Roma group (Ciuchini 
et al., 1993; Ciuchini et ah, 1995). For the other coefficients for which no lattice estimate is available, the following 
"educated guesses" are used: 

• -"3,7,8,9 ~ 

• in the range 1 to 6, in order to account for the large values of b[^1 needed to reproduce the AI = 1/2 rule. 

The parameter Bk is consistently taken from the lattice estimates (Ciuchini et ah, 1995). This determination gives 
in turn the value quoted in Table V for S,, ' by means of the relation Eg ' = Bk which holds if isospin-breaking 

corrections arc ncgicsctccl. 

Finally, because of the matching scale being at 2 GeV, also open charm operators similar to Qi^2 but with the 
strange quark replace by a charm quark (Qf 2) should be included and a value of 2 = -j- 0.15 is assumed. The 
eqs. (4.38)-(4.40) are replaced by 

{Qa)o = {Q3)o + {Q2)o - {Qi)o + {Ql)o - {QDo , (4.41) 

{Q9)o = liQih - liQsh + liQDo , (4.42) 

(<3io)o = {Qi)o + {Q9)o - {Q3)o ■ (4.43) 

The strength of the lattice approach is the direct evaluation of the crucial matrix elements (Qq) and {Qs)2- On 

(2) 

the other hand, while the lattice calculations of B^ ' appear to have settled to reliable numbers, there is still no 
solid prediction for Bq (Gupta, 1998; Martinelli, 1998), and therefore the possibility of sizeable deviations from unity 
remains open. 

The values in table V, which are those used for the current lattice estimate of s'/e , agree with more recent 
determinations (Kilcup ct al., 1998; Gupta ct ah, 1997; Conti ct ah, 1998) except for Bk for which the updated 
central values of 0.92 (Conti et ah, 1998) and 0.90 (Sharpe, 1997) are obtained. 



G. Chiral Quark Model 

Effective quark models of QCD can be derived in the framework of the extended Nambu-Jona-Lasinio (ENJL) 
model of chiral symmetry breaking (For a review, see, e.g.: Bijnens, 1996). Among them is the chiral quark model 
(xQM) (Manohar and Georgi, 1984; Espriu et ah, 1990). This model has a term 

C^QU = -M {qji ^qL + ^ ^Ur) , (4-44) 

added to an effective low-energy QCD lagrangian whose dynamical degrees of freedom are the u,d,s quarks propagating 
in a soft gluon background. The quantity M is interpreted as the constituent quark mass in mesons (current quark 
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masses arc also included in the effective lagrangian). The complete operatorial basis in eq. (2.3) has been analyzed 
for K ^ nn decays, inclusive of chiral loops and complete 0{p^) corrections, by the Trieste group (Bertolini et al., 
1996; Bertolini et al., 19986). 

In the factorization approximation, the matrix elements of the four quark operators are written in terms of better 
known quantities like quark currents and densities, as already shown in eq. (4.3). Such matrix elements (building 
blocks) like the current matrix elements (0| (1 — 75) u \K~^(k)) and (7r"'"(p+)| 57'' (1 — 75) d \K~^{k)) and the matrix 
elements of densities, {0\'s-f5u\K~^{k)) , {Tr~^ {p+)\sd\K'^ (k)) , are evaluated up to O(p^) within the model. The model 
dependence in the color singlet current and density matrix elements appears (via the M parameter) beyond the leading 
order in the momenta expansion, while the leading contributions agree with the well known expressions in terms of 
the meson decay constants and masses. 

Non-factorizable contributions due to soft gluonic corrections are included by using Fierz-transformations and by 
calculating building block matrix elements involving the color matrix T" (see eq. (2.4)): 



I sr 



'T-{l-j,)u\K+{k)) , (7r+(p+)|s7^T"(l-75)d|i^+(fc)) . (4.45) 

Such matrix elements are non-zero for emission of gluons. In contrast to the color singlet matrix elements above, they 
are model dependent starting with the leading order. Taking products of two such matrix elements and using the 
relation 



,a 



9lG%Glp = Y^^GG) {S^J^p - 5,^5,^) (4.46) 

makes it possible to express non-factorizable gluonic corrections in terms of the gluonic vacuum condensate (Pich and 
de Rafael, 1991). The model thus parameterizes all amplitudes in terms of the quantities M, (qq) , and {ugGG /-k) . 
Higher order gluon condensates are omitted. 

The leading order (LO) (0(p°,]3^)) matrix elements {Qi)f'^ and the next-to-leading order (NLO) (0(p^,p^)) cor- 
rections {Qiif^'^ for isospin / = 0, 2 for the pions in the final state are obtained by properly combining the building 
blocks. The total hadronic matrix elements up to O(p^) can then be written: 

{Q,(p))i - [(Q,)j^o + (Q.)f ^°(a*)] + al (m) , (4.47) 

where Qi are the operators in eq. (2.3), and a| (/u) are the contributions from chiral loops (which include wave-function 
renormalization). The scale dependence of the NLO terms comes from the perturbative running of the quark masses. 
The wave-function renormalizations Zk and arise in the xQ^ from direct calculation of the K ^ K and tt — > tt 
propagators. 

The quantities of (/x) represent the scale dependent meson-loop corrections which depend on the chiral quark model 
via the tree level chiral coefficients. They have been included by the Trieste group by consistently applying the MS 
scheme of dimensional regularization. 

At O(p^) the Q5,6 and Qy^g matrix elements contain the NLO coefficients and L%, which within the chiral quark 
model are given by 

^''^-6^fV ("8) 



and 



167r2 24 l&{qq)M \ {qq) 

The hadronic matrix elements are matched with the NLO Wilson coefficients at the scale A,^ ~ 0.8 and the scale 
dependence of the amplitudes is gauged by varying between 0.8 and 1 GeV. In this range the scale dependence of 
e'/e remains always below 15%, thus giving a stable prediction. The 75 scheme dependence, which arise from the 
quark integration in the xQM is also found to numerically cancel to a satisfactory degree that of the NLO Wilson 
coefficients, and the predictions of e'/e in the HV and NDR schemes differ only by 10%. The results reported in the 
following are those of the HV scheme. 

In order to restrict the possible values of the input parameters A'/, {qq) and {agGG/n) , the Trieste group has 
studied the Al = 1/2 selection rule for non-leptonic kaon decay within the xQ^- fitting at the scale /z = 0.8 GeV 
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the calculated amplitudes to the experimental values, they find that within 20% error the A/ =1/2 rule is reproduced 
for 

M = 200 tl MeV , (4.50) 



{asGG/w) = (334 ± 4 MeV)" , (4.51) 

and 

{qq) = (-240t?g MeV)^ . (4.52) 

The fit is obtained for values of the condensates which are in agreement with those found in other approaches, i.e. 

QCD sum rules and lattice, although it is fair to say that the relation between the gluon condensate of QCD sum 
rules and lattice and that of the xQ^ is far from obvious. The value of the constituent quark mass M is in good 
agreement with that found by fitting radiative kaon decays (Bijnens, 1993). 

TABLE VI. The Bi factors in the xQ^- The results for i3i,...,io are shown in the HV scheme, at the scale jj, = 0.8 GeV, for the 
central value of A^^j-, = 340 MeV. The range in the matrix elements of the penguin operators Qs-a arises from the variation 
of (qq) . The value of the (scale and renormalization scheme independent) parameter Bk includes the variation of all input 

paraiiictc^rs. 



Bs 

B4 

Bs ~ Be 

Bf^ ^ Bf^ 

^9 

^(0) 

-"10 

Bf^ ^ Bf^ 
Bk 



9.5 

2.9 

0.41 
-2.3 

1.9 
1.6 ±0.3 
2.5 ±0.1 

3.6 

4.4 
0.92 ± 0.02 

0.41 
1.1 ±0.2 



The obtained factors Bi are given in Table VI in the HV scheme, at /i = 0.8 GeV, for the central value of 
Aqcd (Bertolini et al., 19986). The depende nee on Aqcd enters, as for the Miinchen approach, indirectly via the 
fit of the A7 = 1/2 selection rule and the determination of the parameters of the model. The uncertainty in the 
matrix elements of the penguin operators Qs-s arises from the variation of (qq) . This affects sensibly the Bc^ q 
parameters because of the linear dependence on (qq) of the Qs.e matrix elements in the xQ^) contrasted to the 
quadratic dependence of the corresponding VSA matrix elements. Accordingly, ^5,6 scale as {qq)~^, or via PCAC 
as rriq, and therefore are sensitive to the value chosen for these parameters. It should be remarked that in the xQM 
analysis of (Bertolini et al., 19986) the central value of the quark condensate at the scale n = 0.8 GeV is given 
by (gg)(0.8 GeV) = (-240 MeV)"". As a consequence, the VSA normalization, for the central value of the quark 
condensate, numerically differs from that used in Sect. III.B, which corresponds to {qq){0.8 GeV) = (—222 MeV)'^. 
Finally, it si interesting to notice that decreasing the value of the quark condensate in the xQM depletes the (Qg) 
matrix element relatively to (Qe), and viceversa. 

The parameter Bk is scale and renormalization scheme independent and the error given includes the variation of 
all input parameters (Bertolini et al., 1998a). 

Non-factorizable gluonic corrections are important for the CP-conserving amplitudes (and account for the values 
ofB[°^ and fif^) but are otherwise inessential in the determination of e'/e . 
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H. Discussion 



We would like to make a few comments on the determinations of the matrix elements in the various approaches 
above. 

• All techniques attempt to take into account the AI =1/2 rule, which is the most preeminent feature of the 
physics of hadronic kaon decays. The direct fit of the rule in the phenomenological and lattice approaches 
determines some of the hadronic matrix elements. In the xQ^^ approach, the same fit constrains the few input 
parameters of the model, in terms of which all matrix elements are expressed. The xQ^^ approach is the only 
one for which the fit of the rule determines all hadronic matrix elements. 

Since the operators Qi and Q2, which dominate the AI =1/2 amplitude, do not enter directly the determination 
of e'/e , the way such a fit affects e'/s is indirect and based on the use of operatorial relations as those given in 
eqs. (4.38)-(4.43) in order to obtain information on the matrix elements of some of the penguin operators. 

According to eq. (4.38) a large value of ((52)0 — (Qi)o determines a proportionally large one for ((54)0 if one 
assumes that (Q^jo has a positive value. In the phenomenological approach ((53)0 = 1 is assumed thus obtaining 
a rather large value for B4. Similar values for B4 are obtained, via a similar fit of the A7 =1/2 selection rule, 
in the lattice (see eq. (4.41)). In the xQM, B3 turns out to be large and negative and such that B4 remains 
relatively small, albeit larger than unity. At the same time the value of ((59)0 is increased. The net effect is, by 
looking at the sign of the contributions of the various operators depicted for the VSA in fig. 3, an increase of 
the predicted value of e'/e . 

It would be very interesting to have a lattice estimate of B3 as a test of the two scenarios. 

• The crucial parameters Bq and B^^\ are calculated in the lattice, in the xQ,^^ approach at O(p^), and recently 
by a new estimate of the Dortmund group in 1/Nc at 0{p^) including chiral loops via a cut-off regularization. 

The phenomenological approach varies them according to a 20% uncertainty around their VSA values. 

The xQM finds a substantially larger value for Bq compared to the other approaches. This is due to the meson- 
loop enhancement of the Aq amplitude (Kambor ct al., 1990; Antonelli et al., 1996). It is an open question 
how much of this effect is accounted for in the quenched approximation on the lattice. In addition, the lattice 
calculation of Bq suffers from large renormalization uncertainties. 

The Dortmund group originally found a large enhancement for Bq and suppression for Bg ' . In the latest and 
novel estimate by (Hambye et al., 1998) they find almost no enhancement for B^ and a strong suppression 
for B^K One should wait for a complete 0{p'^) calculation before drawing conclusions from the numerical 
comparison with the xQM results. 

Both the phenomenological approach and the lattice do not include the O(p^) correction terms for the matrix 

elements of the operators Qt.s- The effect of these terms may be within the range of the B^^l values these two 

(2) 

approaclw^s consider. However, when these corrections arc added, they may have the effect of reducing _Bg 
thereby increasing the central value of e'/e . All present calculations of B^ ' agree on a value smaller than the 
VSA result. 

• Those lattice computations that compute the Bi from the K n amplitude, and then obtains the K rnr 
amplitude by means of the chiral lagrangian, use an incomplete 0(jP) lagrangian. In particular, the term 
proportional to G^jj has a vanishing contribution to the K ^ n amplitude, and in order to be determined, the 
knowledge of the K ^ ww amplitude is required. 

• The parameter Bx is numerically the same in the phenomenological and lattice approaches and smaller than 
the xQM result. This parameter has always been a source of disagreement among different estimates. Recent 
lattice determinations tend to assign a larger central value to Bk, closer to the VSA result {Bk = 1)- 

The different values of Bk used in the various approaches lead, as we shall see, to different ranges for the 
relevant combination of CKM matrix elements which enters the determination of e'/e (see section V). 

• The xQ-^-t model approach is the only one for which all matrix elements are actually estimated — and up to the 
O(p^) in the chiral expansion. Of course this approach suffers from its model dependence and the fact that the 
scale and renormalization scheme stability of the computed observables is a numerical feature that is not formally 
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proven (while the lattice and the Miinchen phenomenological estimates are in principle safe in this respect). 
On the other hand, it is the only approach in which the A/ =1/2 rule is well reproduced in terms of natural 
values of the few input parameters when non-factorizable effects like soft-gluon corrections and meson-loops are 

included. These non-factorizable contributions are important in estimating e'/e as shown by the relatively large 
value of Ba and in the interplay between the operators Qi, Q2, Q?, and Q^^ (related by = + Q2 — Qi)- 

• Chiral-loop corrections give in general important contributions to the hadronic matrix elements. A complete 
calculation of the hadronic matrix elements at 0(p*) has been performed only in the framework of the xQM so 
far. 

Of course, it is not sufficient to know the Bi factors in order to predict e'/e , since the impact of the Wilson 
coefficients and other input parameters must also be taken into account. As we shall see, the predictions depend 
crucially on the determination of the relevant CKM entries and the value assigned to rrig (or, via eq. (3.2), the value 
of the quark condensate {qq) ). 



V. THE RELEVANT CKM MATRIX ELEMENTS 

The ratio e'/e , once the measured value of e is used, turns out to be proportional to the combination of CKM 
matrix elements 



Im At = ImVtdVt* , 

which, by using the Wolfenstein parameterization of eq. (1.16), can be written as 

Im Xt:^v>'^A^=r]\Vus\\VcbW 

where A = l^cbl/A^ and A = |Ks|- 

In order to restrict the allowed values of Im At we can solve simultaneously three equations. 

The first equation is derived from eq. (1.13) and gives the constraint from the experimental value of e : 



\Vcbf mS{xt) + mSixx, xt) - r]iS{xc) 



A8 



Bk 



C Alo 



= 0.226 



where 



C 



3V27r2AMLs 



(5.1) 
(5.2) 



(5.3) 



(5.4) 



In writing eq. (5.3) we have neglected in ImA*A( the term proportional to Re At/Re Ac which is of 0{X'^) and used 
the unitary relation ImA* = ImAj. 

Two more equations are those relating 77 and p to measured entries of the CKM matrix: 



1-p 1 



2 I 2 



A^ 
2 



A2 iKbP 
1 \Vtd\^ 

A2 



The allowed values of rj and p are thus obtained, given e , mt, rric and (Barnett et al., 1996) 

\Vus\ = 0.2205 ±0.0018 , 
\Vcb\ = 0.040 ± 0.003 , 
I Kb I /I KbI =0.08 ±0.02. 



(5.5) 
(5.6) 



(5.7) 
(5.8) 
(5.9) 



For \Vtd\ we can use the bounds provided by the measured Bj^-B^ mixing according to the relation (Buras and 
Fleischer, 1997) 



V;w| =8.8-10-3 



200 MeV 

^/Bb^Fb^ 



170 GeV 

mt{mt) 



0.76 



AM 



0.50/ps 



0.5 



0.55 
Vb 



(5.10) 
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The theoretical uncertainty on the hadronic AS = 2 matrix element controls a large part of the uncertainty on 
the determination of ImAt. For the renormalization group invariant parameter Bk we take, as a reference for the 
following discussion, the VSA value with a conservative error of ±30%. 

The AS = 2 parameters r/i 2,3 obtained from QCD are known to the NLO (Buras et al., 1990; Herrlich and Nierste, 
1994; Herrlich and Nierste, 1995; Herrlich and Nierste, 1996). We compute them by taking A^^^^ = 340 ± 40 MeV, 

mb{mb) = 4.4 GeV, mdnic) = 1.4 GeV and mj^°'''^ = 175± 6 GeV, which (in LO) corresponds to mt{mw) = 177 ± 7 
GeV, where running masses are given in the MS scheme. As an example, for central values of the parameters we find 
at fi = rric 

m = 1.33 , 772 = 0.51 , 773 = 0.44 . (5.11) 

This procedure gives two possible ranges for Im At which correspond to having the CKM phase in the I or II quadrant 
{p positive or negative, respectively). Figure 4 gives the results of such an analysis for the central value of mj: the 
area enclosed by the two black circumferences represents the constraint of eq. (5.5), the area between the two gray 
(dashed) circumferences is allowed by the bounds from eq. (5.6); the area enclosed by the two solid parabolic curves 
represents the solution of eq. (5.3) with I3k in the 0.7-^ 1.3 range (notice that the upper parabolic curve corresponds 
to the minimal value of Vcb and vice versa for the lower curve) . 



0.5 
0.4 

0.3 
0.2 
0.1 



-0.4 -0.2 0.2 0.4 

P 

FIG. 4. The allowed rj and p ranges for Bk = 1-0 ± 0.3. 

The gray region within the intersection of the curves is the range actually allowed after the correlation in Vet 
between eq. (5.3) and eq. (5.6) is taken into account. A further correlation is present in going from 77 to Im A^ in 
eq. (5.2). 

In the example of the VSA, where we have taken Bk — 1.0 ± 0.3, from Fig. 1 we obtain 

0.51 X lO"'' < ImAt < 1.6 X 10"'' . (5.12) 

A further refinement of the analysis consists in assigning to each pair of (p, 77) values a Gaussian weight according 
to the deviations from the experimental central values of the computed parameters Vub/Vct, AMb^, s. In this way, a 
Gaussian distribution of the uncertainty on Im At (to be opposed to a fiat one) is found and the error reduced. We 
will use for the discussion of the VSA the flat result of eq. (5.12). 

In general the renormalization group invariant parameter Bk depends on the modeling of the hadronic matrix 
elements, so that different ranges of Im At should be expected according to the different approaches. 

• In the A4iinchen phenomenological approach, where Bk = 0.75 ± 0.15, a range 

0.86 X lO^'^ < ImAf < 1.71 x lO"'' (5.13) 

is found for a flat-distribution of the uncertainties in the input parameters, while the reduced range 

ImAf = (1.29 ±0.22) X 10""* (5.14) 

is obtained for a Gaussian treatment of the same uncertainties. 
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• In the Roma lattice calculation, which takes Bk = 0.75 ± 0.15, the range 

cos (5cp = 0.38 ±0.23 , (5.15) 

is obtained via the Gaussian treatment of the uncertainties, where 5c p is the CKM phase. A result similar to 
eq. (5.14) is found by means of 

|2 \^ub 



ImAt = |K(,|'^Vl-cos2 6cP . (5.16) 

• In the Trieste xQ^ approach, which finds Bk = 1-1 ± 0.2, a flat scan of the input values leads to 

0.62 X lO"'^ < ImAt < 1.4 X lO"'^ . (5.17) 

The larger value of Bk is responsible for the smaller values Im At obtained in this approach. 
For a recent and detailed review on the determination of the CKM parameters see (Parodi et al., 1998). 

VI. THEORETICAL PREDICTIONS 

We have now all the ingredients necessary to understand the various theoretical predictions for e'/e . Let us first 
rewrite eq. (1.11) in such a way that the relationship with the effective operators is more transparent. 
The ratio e'/e can be written as 



— = , , ^ — — Im At 
s 2 e Re^o 



Ho - - n2 



(6.1) 



where, referring to the AS = 1 quark lagrangian of eq. (2.1), 

no = V y^ MQi)o (1 - ^v+v') > (6-2) 

COSdn 

i 

n2 = ^^VyiRe(Qi)2. (6.3) 

cos 02 



(6.4) 



The phase of e'/s is (Maiani et al., 1992) 



77 



-+5o-^2-^. = (0±4)°, (6.5) 

and we can take it as vanishing. Wc assume everywhere that CPT is conserved. An extra phase in addition to (6.5) 
would be present in the case of CPT non-conservation: present experimental bounds constrain it to be at most of the 
order of 10"'' (for a review see (Maiani et al., 1992)). 

Notice the explicit presence of the final-state- interaction phases in eqs. (6.2) and (6.3). Their presence is a 
consequence of writing the absolute values of the amplitudes in term of their dispersive parts. Theoretically, given 
that in eq. (2.2) r <C 1, we obtain 

z, Im(Q,)/ 
J2i Zt Re(Qi)/ 

A phenomenological estimate of the rescattering phases can be extracted from the elastic tt-tt scattering. In chiral 
perturbation theory to 0{p'^) one obtains (Gasser and Meissner, 1991) 

So -621^^2^= 45° ± 6° . (6.7) 

A more recent analysis of pion-nucleon collisions (Chell and Olsson, 1993), based on QCD sum rules and the extracted 
s-wave TT — TT isospin scattering lengths, finds at the kaon mass scale 
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and, accordingly, 



So = 34.2° ± 2.2° 



52 = -6.9° ± 0.2° , 



41° ± 4° . 



(6.8) 



(6.9) 



This result improves on older analyses (Basdevant et al., 1974; Basdevant et al., 1975; Froggatt and Petersen, 1977) 
for which 



So 



37° ± 3° 



-7° ± 1^ 



(6.10) 



All these results are consistent with each other and imply a misalignment of the 7 = over the 7 = 2 amplitude by 

about 20% (cosSq/ cos(52 — 0.8). Final state rcscattcring is not included in the VSA hadronic matrix elements, and in 
the lattice calculations, where the K ^ tt transition is computed. Absorptive components appear when chiral loops 
are included, as in the 1/A^c approach of (Bardeen et al., 1987) and in the xQM approach of the Trieste group. In the 
latter framework the direct determination of the rescattering phases gives at 0{p'^) So — 20° and S2 ~ —12°. Although 
these results show features which are in qualitative agreement with the phases extracted from pion-nucleon scattering, 
the deviation from the experimental data is sizeable, especially in the 7 = component. On the other hand, at O(p^) 
the absorptive parts of the amplitudes are determined only at 0{p^) and disagreement with the measured phases 
should be expected. At any rate, the effect of such a discrepancy on eqs. (6.2)-(6.3) is numerically reduced by the 
cos^o,2 dependence. The authors have therefore chosen to input the experimental values of the rescattering phases in 
all parts of their analysis. This amounts to overstimating systematically the 7 = amplitude by about 15%. Since the 
analysis of the Trieste group is based on the fit of the A7 =1/2 rule with a 20% accuracy, such a bias is reabsorbed 
by the uncertainty found in the determination of (qq). 

Since Im A^^ = according to the standard conventions, the short-distance component of e'/e is determined by the 
Wilson coefficients yi. Because, 2/1 (/x) = 2/2 (m) = 0) the matrix elements of Qi,2 do not directly enter the determination 
of s'/s. 

We can take, as fixed input values: 



2 lei Re An 



349 GeV" 



LO = 1/22.2 . 



(6.11) 



The large value in eq. (6.11) for l/ui comes from the A7 =1/2 selection rule. 

The quantity ilrf+n', included in eq. (6.2) for notational convenience, represents the effect of the isospin-breaking 
mixing between 7r° and the etas, which generates a contribution to A2 proportional to Aq. Clr,+r}' can be written 
as (Donoghue et al., 1986; Buras and Gerard, 1987) 



1 1 



v+v' 



3\/2w 



(cos 0- \/2sm0y + (sin 6* - ^2 cos 6») 



2 ml -ml 



rrir, 



where (Gasser and Leutwyler, 1985) 



mrf - mu 



0.022 ± 0.002 . 



(6.12) 



(6.13) 



The mixing angle 6 has been recently estimated in a model-independent way (Venugopal and Holstein, 1998) to be 

^ = -22° ±3.3°, (6.14) 

which is consistent with the values = —20° ± 4° found in chiral perturbation theory (Donoghue et al., 1986) and 
^ —22° in the 1/A^c expansion (Gasser and Leutwyler, 1985). 
The values above yield 



o n OS +0.03 

— V.ZO _o.o4 ■ 



(6.15) 



Smaller values are found once the uncertainty on the contribution of the rj' is included (Cheng, 1988). For this reason, 
the more conservative range of values used in current estimates of s'/s is 



flri+n' = 0.25 ±0.10 . 



(6.16) 
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A. Toy Models: VSA and VSA+ 



Before summarizing the current estimates of e'/e , it is useful to study some of the steps through which they are 

obtained in a toy model like that given by the VSA. As already pointed out, this model, because of its simplicity, can 
be considered as a convenient reference framework against which all other estimates are compared. 

The VSA+ model that we introduced in Section IV is an attempt to improve on the VSA. It shows how a more 
refined treatment of the electroweak operators, which includes the 0{p^) corrections to the leading constant term, 
can lead to a larger value of e'/e . 

The main purpose of these toy models is to illustrate in a simplified framework some general features of the 
calculation and the impact of some assumptions on the predicted value of e'/e . As we have discussed in Section 
IV, the VSA (as well as the VSA+) cannot give a reliable estimate because of the absence of a consistent scale and 
renormalization scheme matching with the NLO short-distance QCD calculation. 

In the present discussion we use the Wilson coefficients in the HV scheme and set the reference value of the 
matching scale at 1 GeV (see table II). We will then gauge the renormalization scheme dependence of e'/e by varying 
the renormalization scheme from HV to NDR in the VSA amplitudes. Varying the matching scale around 1 GeV will 
show the systematic uncertainty related to the choice of the renormalization scale. 

As we shall sec, different groups work at different renormalization scales because of the peculiarities of their 
approaches. On the other hand, in a consistent approach the choice of the renormalization scale should be immaterial 
as long as observables are concerned. The same holds for the scheme dependence. 

In addition to giving the i?i-parametcrs and the Wilson coefficients in a common scheme and at a common scale, 
one needs to specify the numerical value for the input parameter (qq) which appear in the penguin matrix elements. 
We take the PCAC result, which at 1 GeV and for m„ + = 12 ± 2.5 MeV gives 

{qq) = {-238t{lMeVf . (6.17) 

The mass is often used instead of and m„. Such a change does not reduce the error and may even add further 

uncertainties due to the violations of PCAC that arc larger in the SU{3) case. 

Each of the steps above, necessary in order to estimate e'/e , may carry in practice some model dependence and 
the reader must always bear in mind the assumptions that have entered in the final numerical values. 

Let us now study how the various operators come together to give the final value of e'/e . Figure 5 shows the 
individual contribution of each operator in the VSA (gray histograms) and in the VSA-|- (half-tone histograms). The 
dark histograms show how the various contributions are affected by changing the renormalization scheme from HV 
to NDR in the VSA-I- case. 

The VSA and VSA-|- estimates only differ in the Qr^s matrix elements, as already discussed in section IV. A, while 
moving from HV to NDR affects mostly the Q^^q contributions (see Table II), thus leading to a potentially large effect 
on the VSA prediction for e'/e . 

A central value of the order of 5 x 10"'' is found in the VSA, whereas in going from VSA to VSA-|- the central value 
is increased by 50%. A 25% effect is then related to the renormalization scheme dependence in the VSA+, which 
corresponds to a 50% effect on the VSA. 

Figure 5 shows clearly how systematic effects may sizeably move the e'/e value, due to the change in the destructive 
interference between gluonic and electroweak penguins. 
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FIG. 5. Anatomy of e'/e in the VSA in units of 10 ^ at ^ = 1 GeV with Im At = 1.1 x 10 AU other inputs are taken at their 
central values. Depicted in gray is the VSA, in half-tone the VSA+ and in black the effect of changing the renormalization 
scheme from HV to NDR. 
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FIG. 6. The distribution of the isospin I — (light gray) and 1 = 2 (dark gray) contributions of each operator to e'/e (in units 
of 10"^) in the VSA. 

Figure 6 shows, for the case of the VSA, the distribution of the / = and 2 components in the contributions of 
each operator. This figure is useful in disentangling the role and weight of the individual operators according to the 
isospin projections. 

Finally, in Fig. 7 the value of e'/e in the VSA is shown as we continuously vary the two most relevant parameters: 
Im Xf and (qq) . The two surfaces show in addition the dependence of e'/e on the short-distance input parameters 
Aqcd and rrit as we vary them between their la limits, and include also the dependence on the matching scale which 
is varied from 0.8 to 1.2 GeV. Fig. 7 is useful in showing the correlations between the input parameters and e'/e , 
which qualitatively hold beyond the specific model considered. 
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FIG. 7. Parameter dependence of e'/e in the VSA, in units of 10 . The upper (hghter) surface corresponds to taking /i — 0.8, 
Aqcd ~ 380 MeV and mt = 161 GeV, while the lower (darker) surface corresponds to fj, = 1.2 GeV, Aqcd = 300 MeV and 
mt = 173 GeV. 

From Fig. 7 we can finally extract the range of values taken by the parameter e'/e in the VSA, in the HV scheme, 
as we vary all relevant inputs. Taking into consideration the scale dependence of (qq) we find 

e'/e = (0.5 Igj) X IQ-^ . (6.18) 

Analogously, in the NDR scheme we obtain 

e'/e = (0.8 tJJ) X 10"^ . (6.19) 

The large upper range of the NDR result is a consequence of the increase of the scheme dependence of the Wilson 
coefficients as Aqcd increases and the renormalization scale decreases. 

While the toy models are useful in understanding how various possible contributions enter in the final estimate of 
e'/e , it is clear that some important factors are not included. Among them, the actual range of Im Xt, strictly related 
to the determination of Bk — which might be quite different from the naive VSA — and the consistency of the hadron 
matrix elements with the A/ =1/2 rule — which is important in assessing the confidence level of the e'/e predictions. 
For this reason, we now turn to estimates that incorporate these important features. 

B. Estimates of e'/e 

There are three groups for which an up-to-date calculation is available. In addition we will also briefiy comment 
on some recent partial results obtained within the 1/Nc approach. We will identify the various groups by the names 
of the cities (Miinchen, Roma and Trieste) where most of the group members reside. 

In table VII we collect some of the relevant inputs used by the three up-to-date estimates. There is overall agreement 
on the short-distance input parameters. The Trieste group differs from the other two in the value of Bk, and therefore 
of Im At, that is smaller, and for the inclusion of the FSI phases. The matching scales are different because of the 
peculiarities of each approach which lead to the quoted energy scales. The scale (and renormalization scheme) 
dependence of the final estimates is however rather small. We recall that while this stability is a formal property of 
the lattice and Miinchen phenomenological approaches, it is just a numerical feature of the Trieste estimate. 
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TABLE VII. Comparison of input parameters in various approaches. 



input 


Nliinchen 


Roma 


Trieste 


A (4) 


ozo ± oU MeV 


ooU =h iUU iVlcV 


o4U ± 4U MG V 


mt{mt) 


167 ± D GeV 


157 ± 8 Gov 


Vol ± D GcV 


mb{mi,) 


4.4 (jre\ 


4.0 GeV 


A A C^nAT 

4.4 GeV 


mc{mc) 


l.d GeV 


1.5 GeV 


1.4 GeV 




1.3 GeV 


2 GoV 


0.8 GeV 




150 ± 20 MeV 


128 ± 18 MeV 


220 ± 20 MeV 


{m) (m) 


via PCAC from nis 


via PCAC from rus 


(-240lf° MeV)3 


Bk 


0.75 ±0.15 


0.75 ±0.15 


1.1 ±0.2 


Im At X 10" 


1.29 ±0.22 


1.29 ±0.22 


1.0 ±0.4 


cos (5o 


1 


1 


0.8 


cos52 


1 


1 


1 




0.25 ± 0.05 


0.25 ± 0.10 


0.25 ±0.10 



The experimental value of rtit reported in Table VII — which has become available in the last few years — greatly 
helps in restricting the possible values of e'/e and, as we shall see, rules out, at least for a class of models, any 
mimicking of the superweak scenario by the standard model. 

Starting with eq. (6.1), and given the input parameters in Table VII, the different estimates can be computed by 
writing e'/e in terms of the VSA to the matrix elements and the parameters Be 



2/4-64 + -^2/3-83 



ye Be 



1 



+ 



2V3 



-X 



p 



+ 



■ysBi 



Xj (ysBs + —y7Brj 



1 

n:. 



{ygBg - yioBio) 



(6.20) 



and 



{Qi)2 = ^/6 



p 



a/2 

ysBs + Tr2/7-B7 ) - —X 



N, 



iygBs + yioBio) . 



yrBr + 



(6.21) 



In eqs. (6.20)-(6.21) the values of the parameters L5 and {qq) are obtained according to eq. (4.21) and eq. (3.2) 
respectively, taking into account the scale dependence of the quark masses. 

By inserting the appropriated Bi, taking into account their renormalization scheme dependence, the corresponding 
value of (qq) (or Wg) and the other short-distance inputs, varied within the given uncertainties, the reader can recover 
the estimates for the various groups that are reported in the next few subsections. 



1. Phenomenological Approach 

In the phenomenological approach of the Miinchen group (Buras, Jamin and Lautenbacher, 19936; Buras et al., 

1996) the matching scale is chosen at /i = m,, because it is the scale at which penguins arc decoupled from the CP 

conserving amplitudes and some of the Bi parameters can be extracted from the knowledge of the AI =1/2 rule. 

(2) 

In this approach all Bi except -63,5,6 and B^ are determined from the experimental values of physical processes. 
The operator Q4 rc!c,cives an enhancement due to the rather large value used for B4 that comes from the fit of the 
AI —1/2 rule with the assumption that -63 = 1, as discussed in section IV. E. 

The parameters Bq and B^ ' are varied within a 20% around the VSA values. The quark condensate is written in 
terms of rrig, which is then varied according to the uncertainty of its determination. 
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This procedure yields the two predictions (Buras et al., 1996) 



-1.2 X 10"^ < e'/e < 16.0 x 10"^ , (6.22) 

for TOs(toc) = 150 ± 20 MeV and 

< e'/e < 43.0 x 10"'' , (6.23) 

for ms{mc) = 100 ± 20 MeV. This second range is included in order to study the implications of some recent 
lattice estimates of ms that found such a small values (Gupta and Bhattacharya, 1997; Gough ct al., 1997). Notice 
however that the lower range is somewhat extreme in the light of more recent lattice results now settling down at 
ms{2 GcV) = 110 ± 23 MeV (Bhattacharya and Gupta, 1998) (which corresponds to ms{mc) = 129 ± 27 MeV). 
This range of TOg values is also consistent with recent QCD sum rules estimates (Colangelo et al., 1997; Jamin, 
1998). On the other hand, a substantially larger determination of is obtained from the study of t decays at LEP. 
A preliminary result from the ALEPH collaboration gives ms{mr) = 172 ± 31 MeV (Chen, 1998). It is therefore 
important for the determination of e'/e to understand better the value of this parameter, which via eq. (3.2) affects 
the size of the hadronic matrix elements of the most relevant operators. 

For a Gaussian treatment of the uncertainties that affect the determination of Im Aj, the values (Buras et al., 1996) 

e'/e = (3.6 ± 3.4) x 10"^ , (6.24) 

and 

e'/e = (10.4 ± 8.3) x 10"'' , (6.25) 

are respectively found. 

The same group also gives an approximated analytical formula, in terms of the penguin-box expansion, that is 
useful in discussing the impact in this estimate of the various input values: 

j=lmXtF{xt), (6.26) 

where 

F{xt) = Po + Px Xo{xt) + Py Yo{xt) + Pz Zo{xt) + Pe Eo{xt) . (6.27) 

The a;t-dependent functions in (6.27) are given, with an accuracy of better than 1%, by 

Xo{xt) = 0.660 a;? "^ , Yo{xt) = 0.315 x?'^* , (6.28) 
Zo{xt) = 0.175 xf^^ , Eoixt) = 0.564 a;^-^-^^ . 

The coefRcients Pi are given in terms of Sg^^^' = Bg^^^^mc), Bg^^"^^ = B^^'^\mc) and ms(mc) as follows 



158 MeV 



msirric) + md{mc) 



'(^(6)5(1/2)^^(8)5(3/2)^ . (6.29) 



The Pi must be renormalization scale and scheme independent. They depend however on Aqcd- Table VIII, taken 
from (I 



from (Buras and Fleischer, 1997), gives the numerical values of rf\ rf^ and rf''' for different values of Aqqj-, at 
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TABLE VIII. The penguin-box expansion coefficients for various Aq^^ as given by (Buras and Fleischer, 1997). Only the 
coefficients ro depend at the NLO on the renormalization scheme; the first raw gives their NDR values while the last row shows 
the corresponding values in the HV scheme. The results are given for ms(mc) = 150 MeV. 







A^'cD = 245 MeV 






aW 


= 325 MeV 






A«j3=405 MeV 




i 






i 


i 




rf^ 


(8) 
rl 


i 




i-i 





-2.674 


6.537 


1.111 


-2.747 




8.043 


0.933 


-2.814 


9.929 


0.710 


X 


0.541 


0.011 





0.517 




0.015 





0.498 


0.019 





Y 


0.408 


0.049 





0.383 




0.058 





0.361 


0.068 





Z 


0.178 


-0.009 


-6.468 


0.244 




-0.011 


-7.402 


0.320 


-0.013 


-8.525 


E 


0.197 


-0.790 


0.278 


0.176 




-0.917 


0.335 


0.154 


-1.063 


0.402 





-2.658 


5.818 


0.839 


-2.729 




6.998 


0.639 


-2.795 


8.415 


0.398 



It is important to stress that the approximated formula (6.27), with the numerical coefficient given in Table VIII, 

relics on the values of all Bi used in the phcnomcnological approach, and great attention must be paid to the possible 
effects of the different patterns of Bi and the scale at which they are computed when applying the same formula to 
other frameworks to compare predictions of e'/e in the standard model. 

2. Lattice Approach 

In the lattice approach of the Roma group (Ciuchini et al., 1993; Ciuchini et al., 1995; Ciuchini, 1997), the matching 

scale is taken at /i = 2 GeV. 

As it was for the Miinchen group, the operator Q4 receives an enhancement due to the rather large value used for 
B4 in order to fit the A7 = 1/2 rule with the assumption B3 = 1. The quark condensate is written in terms of mg, 
which is then varied according to the uncertainty of its determination. 

The parameters Bq and Bg ' are explicitly computed on the lattice, although the determination of (Qe) suffers from 
large uncertainties (see section IV. D). 

Only the result obtained via the Gaussian treatment of the errors in the input parameters is reported and yields (Ciu- 
chini, 1997) 

e'/e = (4.6 ± 3.0 ± 0.4) X 10"^ , (6.30) 

where the first error is the variance of the distribution and the second one comes from the residual 75-scheme de- 
pendence. Fig. 8 from (Ciuchini, 1997) shows the anatomy of e'/e in the lattice case. In this figure, the various 
contributions are shown in a manner similar to Fig. 5, with the additional separation of the electroweak components 
in isospin and 2 amplitudes (as in Fig. 6 for the VSA). 
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FIG. 8. Anatomy of e'/e in the lattice approach in terms of the 7 = (A/ = 1/2) and 7 = 2 (A7 = 3/2) components. 

More recent estimates of Bx on the lattice (Kilcup et al., 1998; Gupta et al., 1997; Conti et al., 1998), find a value 
larger than that used in deriving eq. (6.30), which makes Im At and, proportionally, s'/s even smaller. 

3. Chiral Quark Model 

In the xQM approach of the Trieste group (Bertolini et al., 1996; Bertolini et al., 19986), a rather low scale /z = 0.8 
GeV is chosen because of the chiral-loop contribution that become perturbatively too large at scales higher than 
A-^ w nip, the chiral-symmetry breaking scale. Such a low energy scale for the matching makes some of the Wilson 
coefficients larger than in the other approaches and, correspondingly, more sensitive to higher order corrections. 

Let us also recall that the scale and renormalization scheme stability of the computed observable is only a (welcomed) 
numerical feature and no attempt to address formally the cancellation of unphysical dependences is given. On the other 
hand, this estimate is the only one in which all Bi are computed within the same model and in terms of a few basic 
parameters. It is also the only one for which the full 0{p^) amplitudes have been evaluated. It may therefore be quite 
useful in complementing the other estimates by illustrating characteristic patterns of the long-distance contributions. 

The value of Im At is smaller than in the previous two estimates because of the rather large value for Bk (see Table 
III) that is found in this model. 

The quark condensate is a primitive input parameter that is varied according to its determination in fitting the 
Al = 1/2 rule. The value in eq. (4.52) determined at the scale /x = 0.8 GeV by the Trieste group corresponds, via 
PCAC, to Trig — 220 McV. The quark masses appear explicitly in the xQM calculation at the NLO in momentum 
expansion and are treated as independent parameters. It is interesting to observe that in the xQM, because of the 
linear dependence on (qq) of the Qq matrix element, contrasted to the quadratic dependence of {Qs), decreasing 
the value of the quark condensate depletes the destructive interference between the two, and viceversa, partially 
compensating for the overall change of scale. 

Taking into account a 1-a flat distribution of the input parameters the value (Bertolini et al., 19986) 

e'/e = (1.7 11;^) X 10-3 (6.31) 

is found. 

Figure 9 shows explicitly the contributions of the various operators, charted this time operator by operator as in 
Fig 5. 
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FIG. 9. Anatomy of e'/e (in units of 10^'^) within the xQ^ approach. In black the LO results (which includes the non- 
factorizable gluonic corrections), in half-tone the effect of the inclusion of chiral-loop corrections and in light gray the complete 
0(p*) estimate. 

A previous estimate of e'/e by the same group (Bertolini et al., 1996) quoted the smaller value 



The change from (6.32) to (6.31) is due to the following improvements: 

• Inclusion of the complete chiral lagrangian to 0{p^) as discussed in section III. A; 

• Extension of the matrix element calculation to the 0{p^); 

• Update of the short-distance analysis; 

• New ranges of input parameters as determined by the updated fit of the A/ =1/2 rule (Bertolini et al., 1998a). 
4. 1/Nc Approach 

The approach based on a 1/Nc estimate of the hadron matrix elements, including chiral loops, has been first 
pursued by the Miinchen group (Bardeen et al., 1987; Buchalla et al., 1990). Eventually, it was dropped in favor of a 
phenomenological one that was judged to be better. 

Successively, it was taken up by the Dortmund group (Paschos and Wu, 1991; Heinrich et al., 1992; Paschos, 1996). 
Unfortunately, many details of their work are not available and there is no complete and updated calculation. For 
this reason we did not include it in Table VII. 

The latest available estimate quotes the value (Paschos, 1996) 



for nis{l GeV) = 175 MeV. This value is the result of a larger than 1 and a i?g smaller than 1 as obtained by 
including chiral-loop corrections in the matrix elements. 

A very recent and new calculation of Bq and Bg, which addresses systematically the problem of a consistent 
renormalization scale matching between chiral loops and Wilson coefficients, yields a smaller value for Bq while a 

(2) 

much suppressed value for B^ is found (Hambye et al., 1998). No new estimate of e'/e has appeared yet. However, 
some of the relevant observables, as Bk and the 7 = 0,2 amplitudes, show at the present status of the calculation 
a quite poor scale stability (Hambye, 1997; Kohler, 1998), which may frustrate any attempt to produce a reliable 
estimate of e'/e . 



e'/e= (4 ±5) X 10"'* . 



(6.32) 



e'/e = (9.9 ±4.1) x lO"-* , 



(6.33) 
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C. e'/e in the Standard Model: Summary and Outlook 

If we consider that energy scales as different as mt and m^r enter in an essential manner in the determination of 
the ratio e'/e , it is remarkable that this parameter can be predicted at all. Even more remarkable is the fact that all 
theoretical estimates are more or less consistent and a well-defined window of possible values emerges. 

Figure 10 collects the estimates we have discussed and compares them with the two present (1998) experimental 
ranges from CERN (NA31) and FNAL (E731). We have also shown as a reference the results obtained in the simple 
VSA, in the HV and NDR schemes, as discussed at the beginning of this Section. We recall that the VSA error bars 
include a variation of the matching scale from 0.8 to 1.2 GeV. 



Roma Trieste 
(1997) (1997) ; 

Theoretical Predictions 

FIG. 10. Current theoretical predictions for e'/e in units of 10"'^ in the standard model. The horizontal short bars mark the 
central values of each prediction. The two gray areas correspond to the current NA31 and E731 l-cr experimental bounds. In 
dark gray the naive VSA results are shown for comparison (the error bar includes a variation of the matching scale from 0.8 to 
1.2 GeV). 

The two error bars depicted for the Miinchen estimates correspond, from left to right, to flat and Gaussian scanning 
of the input data respectively. Also the reduced size of the error bar of the lattice result is due to the Gaussian 
treatment of the data. 

The entire range between zero and, roughly 3 x 10^"^ is spanned by the available standard model predictions, thus 
dispelling the belief (that has been around in the last few years) that values of the order of 10^"^ were difficult to 
account for within the standard model. 

Given the present theoretical and experimental results it is difficult to draw definite conclusions from their com- 
parison beyond the fact that there are no inconsistencies. On the other hand the forthcoming experimental data may 
crucially help theorists in better understanding the role of non-perturbative QCD in the present estimates. 

To have a pictorial impression of the dramatic improvement expected from the currently running experiments one 
must shrink the experimental ranges within a ± 2 x 10~* error band corresponding to two ticks on the vertical scale 
of Fig. 10. This is shown in Fig. 11 by the horizontal gray band drawn on the central value of the 2a average of the 
NA31 and E731 results 

e'/e = (1.4 ± 1.6) x 10"^ , (6.34) 

which is obtained by following the PDG procedure for error inflation when central values are in disagreement (Winstein 
and Wolfenstein, 1993). 
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FIG. 11. The combined NA31 and E731 experimental bounds (area within the dashed lines), are compared with the most 
updated theoretical estimates for e'/e (in units of 10""^). The gray horizontal band represents the future experimental sensitivity 
shown around the present experimental average value in eq. (6.34). 

Such an improvement in the experimental results will certainly spur a new wave of theoretical analyses. We foresee 
at least three directions along which such a re-analysis may take place: 

• Should the experimental results converge to a common error range of the order of few 10^^, it will be useful 
to focus the attention on the central values obtained by the various theoretical approaches in order to better 
understand the most relevant effects at work. As an example, consider the case that the new experimental 
central value turns out to be near or larger than the present average result. The comparison between the VSA 
and the VSA-I- toy estimates discussed in the present review, together with the results of the Trieste group, 
suggest that the cancellation between the gluon and electroweak penguin operators may be substantially reduced 
once {i) the complete set of the electroweak 0{p^) terms, and (ii) higher order chiral corrections are taken into 
account. These effects can in part be included both in the Miinchen estimate, for those matrix elements that 
are not determined phenomenologically, and in the lattice prediction. Such effects may in fact account for larger 
central values than those presently obtained in those estimates. 

• For what concerns the reduction of the theoretical error, in all estimates a crucial role is played by the knowledge 
of the relevant CKM entries. A large fraction of the theoretical error on e'/e is related to the uncertainty on Im At , 
which amounts by inspection of eqs. (5.13)-(5.16), to a 30%— 40% effect on the total error. The uncertainty 
on ImAt is presently dominated by the determination of Bk- In this respect, a precise determination of ImAt 
from B-physics alone — as expected from the upcoming i3-factories and hadronic facilities — will free this part 
of the analysis from large hadronic uncertainties and thus reduce the impact of non-perturbative QCD in the 
theoretical determination of c'/c . 

• Progress in the lattice estimate of hadronic elements is to be expected in the next few years (Gupta, 1998; Sharpe, 
1998). A reliable estimate of the parameter Be is particularly needed. In addition, achieving the needed precision 
of the order of 10% or below in all relevant matrix elements implies going beyond the quenching approximation. 
The inclusion of higher-order chiral corrections can be also important. Much work is beeing done at present 
which indicates the possibility of addressing quantitatively this issue in the near future. It is from lattice QCD 
that we should expect a conclusive word on the matter. 

VII. NEW PHYSICS AND e'/e 

Physics beyond the standard model may enter the determination of e'/e in many ways. In particular, since the 
origin CP violation is still unclear, 

• It remains an open issue whether the CP violation observed in the K'^-K^ system stems from complex Yukawa 
couplings or from a superweak interaction which goes beyond the standard model; 
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• Even maintaining that the observed CP violation is not superweak in nature, other sources of CP violations 
may be present in addition to, or replace the standard CKM phase in extensions of the standard model. 

• Even if we insist that the CKM phase is the only source of CP violation, new particle contributions to the 
Wilson coefficients of the effective quark Lagrangians still may be relevant for the detailed prediction of CP 
violating observables. 

However, given the discussion of the previous sections and considering in particular the comparison between the 
present theoretical and experimental results shown in Fig. 10, it appears to be a difficult task to disentangle new 
physics effects in e'/e . 

Yet, one question that may be asked is whether the present experimental window allows for visible signals of 
non-standard physics. In order to answer this question we may take the average 2a result of the NA31 and E731 
experiments shown in Fig. 11 and compare it with the range obtained by the union of the most recent theoretical 
estimates (which is a reasonable, albeit biased, procedure). 

It is clear that the case for observable signals of new physics is marginal, to say the least, and that, in order for new 
effects to become visible in e'/e , the next run of experimental data must converge to the most unlikely areas of the 
present range, thus pointing to values of e'/e larger than a few times 10"'^, thereby confirming the 2a upper range of 
the NA31 result, or negative values, thus moving in the lower side of the 2a E731 range. 

For this reason, we think that it is not necessary to present an exhaustive (and exhausting) review of all attempts 
to non-standard physics effects in e'/e . The interested reader may consult (Grimus, 1988; Winstein and Wolfenstein, 
1993; Nir, 1997; Fleischer, 1997) for reviews of possible new-physics effects in CP violation. 

It is nonetheless interesting to analyze whether specific models affect the standard model prediction via definite 
patterns. In order to do so, let us try to infer, inasmuch as possible in a model independent way, how new physics 
may affect the standard model prediction. 

A. Model Independent Analysis 

The key ingredients for a theoretical prediction of e'/e are the determination of Im Af, from the experimental value 
of £ and -B- physics, and the calculation of all direct contributions to e'. These depend, on the short-distance side, 
from the values of the various components of the Wilson coefficients and, on the long-distance side, on the value of 
Bk and the AS = 1 matrix elements for K ttt:. 

If we consider that the new effects modify only the short-distance aspects of the analysis, then the study of e 
exhibits a general feature: the new range of values for Im Xt obtained is always bounded from above by the maximum 
value given in Fig. 4 at p = by the Vub/Vcb measurement, which is a tree level bound and therefore robust to new 
effects. 

As a consequence 

• no sizeable enhancement of e'/e with respect to the standard model estimate can be expected from a modification 
of the short-distance part of e. 

On the other hand, the range of allowed values for Im Xt may be substantially reduced by new physics contributions, 
thus improving on the precision of the e'/e prediction. 

Acting on the matchings of the AS = 1 Wilson coefficients Cj in eq. (2.1) at /z = mw affects the final outcome on 
e'/e . There are patterns on how changing the Ci{mw) may affect the yi at the low energy scale (/x — 1 GeV) via 
strong and electromagnetic renormalizations. 

In Table I, we have schematically reported the distribution of the different types of diagrams that determine the 
initial matching of the Wilson coefficients. Since new heavy particles may show their presence through their virtual 
exchange in the diagrams depicted in Fig. 1, and different type of diagrams show different short distance properties, 
it is important to keep an eye on how the relevant Wilson coefficients are generated 

In Fig. 12 we show examples of how the various coefficients may mix via QCD renormalization and transmit the 
properties of the initial matchings to the other Wilson coefficients at the scale of the low energy process. 
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(a) (b) (c) 




(d.l) (d.2) 



FIG. 12. Effective diagrams showing the one-loop operator mixing via strong renormahzation. 

A direct look at the structure of the LO anomalous dimension matrix of the standard AS = 1 effective quark 
operators is sufficient to show qualitatively how the initial matching conditions may feed down to the final values of 
the various coefficients. 

Here, as a quantitative and model independent test, we have varied the NLO (one-loop) standard-model initial 
matchings Yi = Ci{mw) by factors of (—1, 0, 2) and observed the effects on the corresponding Wilson coefficients yi 
at the scale of 1 GeV. Our conclusions are the following: 

• Only the varying of Y2, Ya, I7 and Yg leads to effects on the low energy j/j larger than a few percents. (Yg and 

matchings remain zero at the one- loop level). 

• Changing the tree level Wilson coefficient Y2, has a proportional effect on all the gluonic penguin coefficients 
(2/3,4,5,6) and similarly on yii,i2, because of the large additive renormahzation induced via the insertion of Q2 
in the penguin like diagrams (d) in Fig. 12. The influence on yg of changing Y2 by a few ten percents is 
therefore dramatic for the prediction of e'/e . On the other hand, one needs a new particle replacing tree level 
W exchanges and tree level physics constrains dramatically these contributions. It is therefore unlikely to expect 
sizeable deviations of Y2 from its standard model value. 

• Changing Ig itself in the range given has no much effect on which is affected always less than 10%, and it 
affects j/8 at the percent level. Multiplicative renormahzation is not the leading renormahzation effect for the 
gluonic penguins. 

• Changing Y7 modifies proportionally t/7 and and may have therefore a dramatic impact on e'/e . 

• Changing Yg modifies proportionally yg and yio and may affects e'/e at the few 10% level via the contribution 

of Qg. 

It seems therefore that the most relevant potential for new physics effects on e'/e resides in the electroweak penguin 
sector (see Table I). As a matter of fact (Buras and Silvestrini, 1998) have recently shown that bounding the 
contribution of the effective sdZ vertex via e'/e leads to the strongest constraints on some rare kaon decays which are 
governed by Z— penguin diagrams. 

On the other hand, new physics modifications of the standard-model penguin and box diagrams for AS' =1,2 
transitions affect also the corresponding AB = 1,2 amplitudes. It is therefore likely that in a specific model the 
experimental bounds coming from B physics may indirectly constrain the deviations on the electroweak initial match- 
ings within a few 10% (Nir, 1997; Fleischer, 1997). These bounds would make it hard for new physics to show up in 
visible deviations from the standard e'/e prediction. 

The past literature on the subject confirms the general conclusion that we reached in the above discussion. The 
effect on e'/e of charged Higgs particles in the two Higgs model has been studied (Buchalla ct al., 1991). The same 
problem has also been discussed in the more general framework of softly broken supersymmetry (Gabrielli and Giudice, 
1995). In both cases no significant departures from the standard model are expected once all bounds are properly 
implemented. 
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VIII. CONCLUSIONS (MARCH 1999) 

On February 24, 1999 the KTeV collaboration has announced ^ a preliminary result based on the analysis of 20% 
of the data collected which gives 



Re e'/e = (28 ± 3.0 (stat) ± 2.6 (syst) ± 1.0 (MC stat)) x 10" 



(8.1) 



This result largely deviates from the previous E731 value of eq. (1.20) and sits in the ballpark of the NA31 result (1.19). 
This value of e'/e , if confirmed, signals with high confidence the presence of direct CP violation in kaon decays, closing 
successfully a longstanding and challenging experimental quest. Theoretically, the superweak scenario (Wolfenstein, 
1964) is then excluded as the sole source of CP violation. 

Averaging (8.1) with (1.19) and (1.20), together with the older E731 result Re e'/e = (32 ± 30) x 10""^, leads to 
the value 



Re e'/e = (21.8 ±3.0) x 10" 



■2) 



In Fig. 13 we update the comparison between theory and experiment including the new KTeV result. The light 
gray area shows the 2a range of the average in eq. (8.2). 



e' /ex 10^ 



KTeV 



NA31 +E731 + KTeV 




SNA31 +E731 



Muenchen Roma Trieste 



Experiments vs. Theory 

FIG. 13. The new 2a KTeV result (area enclosed by the long-dashed lines) is compared with the combined 2a experimental 
bounds of NA31 and E731 (area enclosed by the short-dashed lines). The combined average in eq. (8.2) is shown {2a) by the gray 
band. The Miinchen, Roma and Trieste theoretical estimates for e'/e are shown with their central values. The Miinchen and 
Roma predictions include gaussian treatment of the input parameters while the uncertainty in the Trieste estimate corresponds 
to flat parameter spanning. The second Miinchen prediction (light gray) corresponds to taking a low rris range (see eq. (6.25)). 

The comparison between the present experimental average and the theoretical predictions shows a substantial 
deviation from the Roma and Miinchen estimates. Such a disagreement may be reduced by considering a light rrig 
(see the discussion after eq. (6.22)). The two Miinchen predictions shown in Fig. 13 correspond to mg{mc) — 150 ±20 
and ms{mc) = 100 ± 20 (light gray), respectively. 

The Trieste estimate is rather insensitive to TOs since this parameter enters explicitely only at the NLO in the chiral 
expansion, while the value of the quark condensate is determined by the fit of the A/ =1/2 selection rule. 

Let us recall that the most recent (quenched) lattice estimates of nis find ms(2 GeV) = 110±23 MeV (Bhattacharya 
and Gupta, 1998), corresponding to ms{mc) — 129 ± 27 MeV. This range of is also consistent with recent QCD 
sum rules estimates (Colangelo et al., 1997; Jamin, 1998), while a substantially larger determination of iris is obtained 
from T decays at LEP. A preliminary result from the ALEPH collaboration gives nis{mr) = 172 ± 31 MeV (Chen, 
1998). In order to assess the theoretical implications of the KTeV result it will be important to better understand 



^See http : //f nphyx-www. fnal . gov/experiments/ktev/ktev .html 
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the value of nis , that is currently used to parametrizee the hadronic matrix elements of the crucial operators Qe and 
Qs (for a more detailed discussion see (Keum et al., 1999)). 

On the other hand, as we argue in the summary of section VI, it is premature to take a value of s'/s at the 10~^ 
level as a signal of new physics. In particular, it is worth observing that: 

• One of the standard model predictions, which via the A/ = 1/2 selection obtains all matrix elements, is in good 
agreement with the experimental average in eq. (8.2); 

• As it is shown by the VSA and VSA+ toy models, and as it appears from the Trieste calculation, the inclusion 
of NLO order chiral corrections might alone conspire toward increasing the standard model value obtained in 
the present phenomenological and quenched lattice predictions. 

At any rate, efforts on improving all theoretical estimates are now required. In particular, a confident assessement 
of the size of Bq from lattice will be of crucial relevance. 

As discussed in section VI. C, the uncertainty in all present theoretical estimates may be substantially reduced by a 
better determination of ImAt, whose error is presently dominated by the uncertainty on Bk- A precise determination 
of Im At is expected in the upcoming years from i?-physics at the B-factories and at the hadronic colliders. In addition, 
the rare kaon decays 'k'^vv and —> n^vD provide together a clean probe of ImA* (Buchalla and Buras, 

1996). While the latter may be seen at the Brokhaven National Laboratory within the next year, a new experiment 
at the same laboratory has been approved to measure Br(i^L w^vD) with a 10% precision by the year 2005. This 
will allow a determination of Im Af with a similar accuracy. 

In conclusion, the determination of e'/e is a great challenge to both experimentalists and theorists. As more precise 
experimental data become available, improvements in the theoretical calculations are also expected. The interplay of 
the two will hopefully shed more light on the flavor structure of the standard model and on some non-perturbative 
aspects of QCD. 
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